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Ao se falar em Educação a Distância (EAD), é preciso ressaltar a possi- 
bilidade que essa modalidade de ensino promove de inserção social 
pela disseminação de meios e processos de democratização do conheci- 
mento. A meta é não só elevar os índices de escolaridade, mas também 
oferecer uma educação de qualidade, disponibilizando uma formação 
inicial e/ou continuada, em particular, a professores que não tiveram 
acesso a esse ensino. 


Não se pode ignorar que é fundamental haver, sempre, plena conexão 
entre educação e aprendizagem. A modalidade a distância é um tipo 
de aprendizagem que, em especial na Universidade Federal de Minas 
Gerais (UFMG), já está concretizada como um ensino de qualidade. 
Hoje, a aprendizagem tornou-se, para todos os profissionais dessa 
universidade envolvidos no programa de Educação a Distância, sinô- 
nimo de esforço e dedicação de cada um. 


Este livro objetiva desenvolver conhecimentos essenciais para o seu 
estudo. Os alunos estudarão o material nele contido e muitos outros, 
que lhe serão sugeridos em bibliografia complementar. É importante 
terem em vista que essas leituras são de extrema importância para, 
com muita dedicação, avançarem em seus estudos. 


Cada volume da coletânea está dividido em aulas e, em cada uma delas, 
trata-se de determinado tema, que é explorado em diferentes formas 
— textos, apresentações, reflexões e indagações teóricas, experimenta- 
ções ou orientações para atividades a serem realizadas pelos alunos. Os 
objetivos propostos em cada uma das aulas indicam as competências e 
habilidades que os alunos, ao final da disciplina, devem ter adquirido. 


Espera-se que, assim, eles se tornem autônomos, responsáveis, críticos 
e decisivos, capazes, sobretudo, de desenvolver a própria capacidade 
intelectual. Os alunos não podem se esquecer de que toda a equipe 
de professores e tutores responsáveis pelo curso estará, a distância ou 
presente nos polos, pronta a ajudá-los. Além disso, o estudo em grupo, 
a discussão e a troca de conhecimentos com os colegas serão, nessa 
modalidade de ensino, de grande importância ao longo do curso. 


Agradeço aos autores e à equipe de produção pela competência, pelo 
empenho e pelo tempo dedicados à preparação deste e dos demais 
livros do EAD. Espero que cada um deles possa ser valioso para 
os alunos, pois tenho a certeza de que vão contribuir muito para o 
sucesso profissional de todos eles, em seus respectivos cursos, na área 
da educação em geral, no País. 


Ione Maria Ferreira de Oliveira 

Coordenadora do Sistema Universidade Aberta do Brasil (UAB/UFMG) 
Assessora do Centro de Apoio de Educação a Distância (CAED) 

de setembro de 2009 a abril de 2010 
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Apresentação 


Este é um livro para uma disciplina introdutória de Equações Diferen- 
ciais Ordinárias destinado a alunos do curso de licenciatura em Mate- 
mática a distância. 


Ele é dividido em seis aulas. No início da Aula 1 é feita uma introdução 
às equações diferenciais em geral. Depois, entre as equações de 1? 
ordem, são estudadas as equações lineares. Na Aula 2 são trabalhadas 
as equações separáveis. Terminamos a aula com algumas aplicações. 
Na Aula 3 são vistas a análise qualitativa de equações autônomas e a 
existência e unicidade de soluções de equações diferenciais. 


As equações lineares de 2º ordem é o assunto das aulas 4, 5 e 6. Nas 
aulas 4 e 5 o estudo, tanto das equações homogêneas como das equa- 
ções não homogêneas, é feito inicialmente no caso geral, e depois 
no caso particular, em que os coeficientes são constantes. Na Aula 6 
temos como aplicações o estudo das oscilações. 


Os desenhos e gráficos deste livro foram feitos usando o MATLABº* 
com o pacote GAAL, disponível na web no site do autor e o Maxima 
também com o pacote GAAL também disponível no site www.mat. 
ufmg.br/-regi do autor desta obra. Nela também estão disponíveis 
páginas interativas. 


Gostaria de agradecer às professoras Joana Darc A. S. da Cruz e Maria 
Cristina C. Ferreira pela revisão do presente texto e pelas sugestões 
apresentadas. 


O autor 


*MATLAB é marca registrada de The Mathworks, Inc. 


Exponenciais e Logaritmos 


x = e™* = Inte. 
etb = e'e eln(ab) = lna +1nb. 


In £ = Ina — Inbe blna = Ina. 


Números Complexos 


Pré-requisitos 


Se z = a + ib € C, então i? = —1,a parte real e a parte imaginária de z são os números reais 


R{z} =a é S{z} =b. 


Se z4 =a + ib,z2 = c + id € C, então 


zı +z2=(a+0)+i(b+d), zız2 = (ac — bd) + i(ad + bc). 


As raízes da equação do 2º grau ax? + bx + c = 0 com A = b? — 4ac < O são 


x = = 
oa 2a | 


Trigonometria 


cos(a + b) = cos a cos b F sen asen b. 


sen (a + b) = sena cos b + sen b cosa. 


Derivação 


d 
sen x = cos x, y COSxX = —sen x. 
x 


dx 


d u(x) — du u(x) 
Fra = Fria $ 
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Integração 


° fyo) + Bz(x))dx = a | y(x)dx +B [alada 
[suma = [reay. 


° fudo=uo- f vdu. 


n+1 
° fax a +c. 


n+1 


° fsenxar = -cosx +c, | cosxdx — sen +c. 


AX 


e 
/ e™dx = — +c. 
a 
Decomposição em Frações Parciais 


p(x) A Ak 


= e . Aqui É k. 
EE E = qui o grau de p(x) é menor do que 
A A B B 

p(x) EE. Sed m p l a "2 —, Aqui o grau de 
a(x — a1)" (x — ag) x-a (x— a)" x-a (x — &2)"2 
p(x) é menor do que n + nz. 

p(x) SA eua An . Bco E Ba, + CX 

a(x — a1)" (x2 + ax Bo) x-a (x= a)" x2 + ax + b2 (x2 + ax +B: 


Aqui o grau de p(x) é menor do que nı + 2n2 e A = a4 — 48, < 0. 
Matrizes e Sistemas Lineares 


Uma matriz quadrada é invertível se, e somente se, o seu determinante é diferente de zero. 


Um sistema de equações lineares 


ax + by = e 
cx + dy = 


l 
s“ 


pode ser escrito na forma 


a b x e 
aG pe 
Se a matriz do sistema A é invertível, então para todo par (e, f) o sistema tem uma única solução 
(x,y) (A solução é X = A-1B). 


em que 


Equações Diferenciais de 1º Ordem 
1º Parte 


1.1 Introdução às Equações Diferenciais 


OBJETIVOS 
Ao terminar esta seção você deverá ser capaz de: 


e Compreender o que é uma equação diferencial e o que é solução de uma 
equação diferencial. 


e Classificar uma equação diferencial quanto ao tipo, a ordem e a linearidade. 


e Saber o que é uma solução de uma equação diferencial de 1º ordem e de um 
problema de valor inicial. 


Uma equação algébrica é uma equação em que as incógnitas são números, 
enquanto uma equação diferencial é uma equação em que as incógnitas 
são funções e a equação envolve derivadas destas funções. Numa equação 
diferencial em que a incógnita é uma função y(t), t é a variável indepen- 
dente e y é a variável dependente. Vejamos alguns exemplos. 


Exemplo 1.1. O movimento de um pêndulo simples de massa m e compri- 
mento 1 é descrito pela equação diferencial 

do g 

JE + T sen 0 = 0. 
Nesta equação a incógnita é a função 0 (t). Assim, 0 é a variável dependente 
e t é a variável independente. 
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Figura 1.1 - Pêndulo Simples 


F =-yv F xt Focos(mt) 
-— 
F=-kx 
e 
F==yv Fot = F coslot) 
F=-yv 


Fg = F cos(ot) 
+" TD" +» 
F=-kx 
e 


x 


Figura 1.2 - Sistema massa-mola 


Exemplo 1.2. Um sistema massa-mola composto de um corpo de massa 
m preso a uma mola com constante elástica k, sujeita a uma força de re- 
sistência F, = —yv = —y% e uma força externa Fex(t) = Focos(wt), é 
descrito pela equação diferencial 
a dx à dx 
de Va 


+kx = Focos(wt). 


Nesta equação a incógnita é a função x(t). Assim, x é a variável depen- 
dente e t é a variável independente. 


Exemplo 1.3. Numa região do plano em que não há cargas elétricas o 
potencial elétrico u(x, y) em cada ponto (x,y) da região satisfaz a equação 
diferencial 

du 


du 
ox? 


dz é 


AULA 1 


Nesta equação a incógnita é a função u(x,y). Assim, u é a variável depen- 
dente, e x e y são as variáveis independentes. 


R 


V(t) 
Figura 1.3 - Circuito RC 


Exemplo 1.4. Um circuito RC é um circuito que tem um resistor de 
resistência R, um capacitor de capacitância C e um gerador que gera uma 
diferença de potencial V(t) ligados em série. A carga Q(t) no capacitor é 
descrita pela equação diferencial 

dQ 1 

R= + =Q = V(t). 

E oger 
Nesta equação a incógnita é a função Q(t). Assim, Q é a variável depen- 
dente, e t é a variável independente. 


1.1.1 Classificação 


As equações são classificadas quanto ao tipo, a ordem e a linearidade. 


1. Quanto ao tipo uma equação diferencial pode ser ordinária ou 
parcial. Ela é ordinária se as funções incógnitas forem funções de 
somente uma variável. Caso contrário, ela é parcial. Portanto, as 
derivadas que aparecem na equação são derivadas totais. Por exem- 
plo, as equações que podem ser escritas na forma 


F(t,y,y',y",..) =0, 


em que y é função apenas de t, são equações diferenciais ordinárias, 
como as equações dos Exemplos 1.1, 1.2 e 1.4. A equação do Exemplo 
1.3 é parcial. 


2. Quanto à ordem uma equação diferencial pode ser de 1°, de 2º, ..., de 
n-ésima ordem, dependendo da derivada de maior ordem presente 
na equação. Uma equação diferencial ordinária de ordem n é uma 
equação que pode ser escrita na forma 


F(t,y, y, y", y ™) =0. 


As equações dos Exemplos 1.1,1.2 e 1.3 são de 2º ordem, e a equação 
do Exemplo 1.4 é de 1º ordem. 
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3. Quanto à linearidade uma equação diferencial pode ser linear ou 
não linear. Ela é linear se as incógnitas e suas derivadas aparecem 
de forma linear na equação, isto é, as incógnitas e suas derivadas 
aparecem em uma soma em que cada parcela é um produto de al- 
guma derivada das incógnitas com uma função que não depende das 
incógnitas. Por exemplo, uma equação diferencial ordinária linear de 
ordem n é uma equação que pode ser escrita como 


ao(t)y + TOE Haz(t) m reid an(t) Jr = f(t). 


As equações diferenciais ordinárias que não podem ser colocadas 
nessa forma são não lineares. As equações dos Exemplos 1.2, 1.3 e 
1.4 são lineares, e a equação do Exemplo 1.1 é não linear. 

1.1.2 Soluções de Equações Ordinárias 

Uma solução (particular) de uma equação diferencial ordinária de ordem 

n em um intervalo I é uma função y(t) definida no intervalo I tal que as 


suas derivadas de ordem até n estão definidas no intervalo I e satisfazem 
a equação neste intervalo. 


Exemplo 1.5. Considere a equação 
ay" +by +cy=0, coma,b,c € R,a £ O tais que b? — 4ac = 0. 


Vamos mostrar que y(t) = em! é solução desta equação. 


Substituindo-se y(t), y'(t) e y”(t) no primeiro membro da equação, 
obtemos 


b? b 
ay" +by +cy = age b( e-h) 4 ce~ zt 
o b? o b? 4 ejeêt 
4a 2a 
-b +4 ) 
+4ac bro 0, 
4a 


pois, por hipótese, b? — 4ac = 0. Assim, y(t) = e =! é solução da equação. 


1.1.3 Equações Ordinárias de 1º Ordem 


As equações diferenciais ordinárias de 1º ordem são equações que podem 
ser escritas como 


F(t,y,y) =0. 
Vamos estudar equações de 1º ordem, que podem ser escritas na forma 


= (ty) (11) 


AULA 1 


Uma solução (particular) de uma equação diferencial (1.1) em um inter- 
valo I é uma função y(t) definida no intervalo I tal que a sua derivada y’ (t) 
está definida no intervalo I e satisfaz a equação (1.1) neste intervalo. 


O problema 
dy 
I fit 
ar Te (1.2) 
y(to) = yo. 
é chamado problema de valor inicial (PVI). Uma solução do problema de 
valor inicial (1.2) em um intervalo 1 é uma função y(t) que está definida 


neste intervalo, tal que a sua derivada também está definida neste intervalo 
e satisfaz (1.2). 


Quando resolvemos uma equação diferencial ordinária de 1º ordem obte- 
mos uma família de soluções que dependem de uma constante arbitrária. 
Se toda solução particular puder ser obtida da família de soluções que 
encontramos por uma escolha apropriada da constante, dizemos que a 
família de soluções é a solução geral da equação. 


Exemplo 1.6. A equação 
dy _ 3t 
P 

pode ser resolvida por integração direta obtendo 


- fam 
vt)= [e dt = 5 +€, 


que é a solução geral da equação diferencial dada. 


Para encontrarmos a solução do PVI 


dy 3 
e 
v(1/3) = e/3 


substituímos t = 1/3 e y = e/3 na solução geral encontrada, obtendo 
C = 0. Assim, a solução do PVI é 


y(t) = 3" 


válida para —oo < t < œ, que é o maior intervalo em que a solução e sua 
derivada estão definidas. 
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Figura 1.4 - Soluções da equação e do PVI do Exemplo 1.6 
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AULA 1 


. Classifique as equações abaixo quanto ao tipo, a ordem e a linearidade. 
(a) yy' +t=0. (b) x?y" + bxy’ +cy =0. 


. Determine qual ou quais das funções yı (x) = x2, y2(x) = x? e y3(x) = e-* são soluções da 
equação 


. Sejam a,b,c € R. Mostre que 


(a) y(t) = e, com r raiz de ar + b = 0, é solução da equação ay' + by = 0. 
(b) y(t) = e", com r raiz de ar? + br + c = 0, é solução da equação ay” + by' + cy = 0. 
(c) y(x) = x", comr raiz de r? + (b — 1)r +c = 0, é solução da equação x?y” + bxy' + cy = 0. 


. Determine os valores de r para os quais a função y(t) é solução das equações: 

(a) y(t) = — ey +ty? =0. OMO z ey —6ty = 0. 
LR ya Did = ses 
b) y) = z gey -2y =o. (O y() = 5 


ey — ty? =0. 
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1.2 Equações Lineares de 1º Ordem 


OBJETIVOS 
Ao terminar esta seção você deverá ser capaz de: 


e Identificar uma equação diferencial linear de 1º ordem. 
e Calcular o fator integrante de uma equação diferencial linear de 1º ordem. 
e Encontrar a solução geral de uma equação diferencial linear de 1º ordem. 


e Resolver um problema de valor inicial correspondente a uma equação diferencial 
linear de 1º ordem. 


As equações (diferenciais ordinárias) lineares de 1º ordem são equações 
que podem ser escritas como 


d 

Z + p(t)y = q(t). (1.3) 
1.2.1 Equações em que p(t) = 0 

Se a função p(t) = 0, a equação (1.3) torna-se 


d 
Z = q(t), (1.4) 


que é facilmente resolvida integrando-se os dois lados. Assim, a solução 
geral desta equação é dada por 


y(t) = fatnát+c. 


Exemplo 1.7. A equação 
d 
= = sen (2t) 
pode ser resolvida por integração direta obtendo-se a solução geral 


yŒ) = f sen(2t) dt = =. a 


Na subseção 1.2.2 e na seção 1.3 veremos técnicas de se encontrar soluções 
de equações de 1? ordem que se baseiam em transformar a equação inicial 
em uma equação do tipo (1.4). 


1.2.2 Equações Lineares - Caso Geral 


Vamos considerar equações da forma 


d 
a T POY = a(t). (1.5) 


AULA 1 


Figura 1.5 - Soluções da equação do Exemplo 1.7 


Vamos definir uma função auxiliar, u(t), de forma que, ao multiplicarmos 
a equação (1.5) por esta função, a equação obtida seja uma equação linear 
com p(t) = 0, ou seja, do tipo (1.4), que já resolvemos anteriormente. 
Uma função com esta propriedade é chamada fator integrante da equação 
linear. 


Seja 
u(t) = ef pltdt, 


Vamos mostrar agora que u(t) = ef P(Ddt é um fator integrante da equação 
(1.5). 


Observe em primeiro lugar que 


= elrone ( f pdt) PO) =p (19) 


Assim, multiplicando-se (1.5) por u(t), obtemos 
d 
uE + (bp ly = uO. (1.7) 


Mas, como por (1.6), u(t)p(t) = a então (1.7) pode ser reescrita como 


d d 
nE + Ey = pa). (1.8) 


Entretanto, o lado esquerdo desta equação é a derivada de um produto, o 
que faz com que ela possa ser reescrita na forma 


E WOYE) = uA (19) 


A equação (1.9) é uma equação do tipo (1.4), ou seja, 
dY 
TO 


em que Y(t) = u(t)y(t) e f(t) = u(t)q(t). Assim, a solução geral de (1.9) é 
dada por 


uO = [ula(nar+c. 
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Como u(t) £ 0, para todo t € R, dividindo-se a equação anterior por u(t), 
obtemos que a solução geral de (1.5) é dada por 


yO = zp ( f OOd +c). 


Mostraremos na subseção 1.2.3 como podemos chegar a u(t) = e! ptdt 
como fator integrante da equação (1.5). 


Atenção: Não se deve memorizar a fórmula obtida no final. O que fize- 
mos aqui foi mostrar o caminho que deve ser seguido para resolver uma 
equação linear de 1? ordem. 


No próximo exemplo vamos seguir os mesmos passos que seguimos no 
caso geral. 


Exemplo 1.8. Considere a equação 


O fator integrante é 


u(t) eJ 2dt ent ein? £. 


Multiplicando-se a equação acima por u(t), obtemos 


O lado esquerdo é igual à derivada do produto t2y(t). Logo, a equação 
acima é equivalente a 


df, 3 
a (Pye) =t 
Integrando-se, obtemos 
A 
2 == Es 
t y(t) =7 +C 


Explicitando y(t), temos que a solução geral da equação diferencial é 


PC 
y(t) = ITE (1.10) 


Podemos esboçar as soluções desta equação diferencial. Para C = 0, a 
solução é a parábola 
2 
H=. 
= 


Para C Æ 0, temos que o domínio de y(t) é o conjunto dos números reais 
tais que t Æ 0. limp-,+oo y(t) = +00, se C £ 0. Além disso, 


limy(t) = +œ, seC>0 


t>50 


limy(t) = —oo, seC <0. 
t50 


AULA 1 


Vamos analisar o crescimento e decrescimento das soluções 


dy t 20 q 
d 2 Po” 
se, e somente se, 
tt =4C. 


Assim, se C > 0, as soluções têm somente pontos críticos em t = +V4C e 
se C < 0, elas não têm ponto crítico. 


Exemplo 1.9. Considere o problema de valor inicial 


dy 2 
ar Fay! 


y(2) =3. 


A equação é a mesma do Exemplo 1.8. Substituindo-se t = 2 e y = 3 em 
(1.10), obtemos 
J= 4 n C 
4 4 
De onde obtemos que C = 8. Portanto, a solução do problema de valor 
inicial é à 
= 8 
=>. 
Observe que a solução deste problema de valor inicial é válida no intervalo 
(0, +00), que é o maior intervalo contendo t = 2 (pois a condição inicial é 
y(2) = 3), em que a solução e sua derivada estão definidas. Se a condição 
inicial, ao invés de y(2) = 3, fosse y(—2) = 3, a solução teria a mesma 
expressão, mas o intervalo de validade da solução seria (—o0,0). 


Figura 1.6 - Soluções da equação do Exemplo 1.7 
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1.2.3 Como chegar ao fator integrante u(t) = e! POd ? 


Vamos mostrar como podemos chegar ao fator integrante p(t) = e! Pdt, 
Comparando-se as equações (1.7) e (1.8) na página 21 vemos que o fator 
integrante u(t) deve ser uma função que satisfaz a equação diferencial 


= p(Dy(t) 


Esta é também uma equação linear, mas com q(t) = 0. Supondo-se u(t) £ 
0, vamos multiplicar esta equação por 1/ u(t), obtendo a equação 


1 du 
uE de p(t). 
Como zo = Fm (In | (t)|), a equação anterior pode ser reescrita como 


a ON E = po. 


Mas pela Regra da Cadeia esta equação é equivalente a 
z mO) = 
AET P 


que é uma equação do tipo (1.4) que pode ser resolvida simplesmente 
integrando-se ambos os membros obtendo 


in[u()| = f pd +C. 


Aplicando-se a exponencial a ambos os membros e eliminando-se o valor 
absoluto, obtemos 


u(t) = erel PlDdt — Cel plDat, 


Como estamos interessados em apenas um fator integrante, podemos to- 


mar C = 1, e obtermos 
u(t) zZ ef pltjdt, 


AULA 1 


Exercícios 


. Resolva os problemas de valor inicial: 


'+(1-2x)y = xe* y’ — costy = tef tsent 
z (9 E 
y(0) = y(0) = 
-B 5 
b) vi ! tay- B+ (d) y' +xty= xte 5 
y(0) = (0) = 1. 
Resolva as equações 
A (c) y — fy = we 
6) y= 
. (a) Resolva o problema de valor inicial: 
À y E 5xty = xó 
y(0) = yo. 
(b) Para quais valores de yọ a solução é crescente e para quais valores de yọ a solução é 


decrescente. 
(c) Qual limite de y(x) quando x tende a +0. O limite depende de yo? 


. (a) Resolva o problema de valor inicial: 
À (x? = 0 +xy =0 
y(5) = yo. 
(b) Qual o intervalo de validade da solução? 
(c) Qual limite de y(x) quando x tende a +o. O limite depende de yọ? 


. Considere a equação 


dy = 
mT p(t)y =0. 
(a) Mostre que, se y1 (t) e y2 (t) são soluções da equação, então y(t) = y1 (t) + y2 (t) também 


oé. 
(b) Mostre que, se y1 (t) é solução da equação, então y(t) = cy (t) também o é, para qual- 
quer constante c. 


. Considere as equações 


dy = 

T +p(t)y =0, (1.11) 
dy = 
ar POY = q(t). (1.12) 


Mostre que, se y1(t) é solução da equação (1.11) e y>(t) é solução da equação (1.12), então 
y(t) = cy1 (t) + y>(t) é solução de (1.12), para qualquer constante c. 
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Equações Diferenciais de 1º Ordem 
2° Parte 


2.1 Equações Separáveis 


OBJETIVOS 
Ao terminar esta seção você deverá ser capaz de: 


e Identificar uma equação diferencial separável. 
e Encontrar a solução geral de uma equação diferencial separável. 


e Resolver um problema de valor inicial correspondente a uma equação diferencial 
separável. 


As equações (diferenciais ordinárias) separáveis são equações que podem 
ser escritas na forma 


d 
sl) Te = SO) (2.1) 
Seja 


então 


dh 
dy ~ 80). 
Eca dh : 
Substituindo-se g(y) por Er na equação (2.1), obtemos 
dh dy 
dy dx = f(x). (2.2) 


Mas, pela Regra da Cadeia, 


d — dhdy 
moO = Gado 


o que implica que (2.2) pode ser escrita como 


d 
EM(U()) = fla). 23) 
A equação (2.3) é do tipo (1.4), página 20, ou seja, é da forma 
d 


Y 
dx = f(x), 
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em que Y(x) = h(y(x)). Assim, integrando-se (2.3) dos dois lados, obte- 
mos que a solução geral de (2.1) é dada implicitamente por 


nya) = f Faar+e. 


Também podemos obter a solução da maneira mostrada a seguir. Integrando- 


-se em relação a x ambos os membros de (2.1), obtemos 


feo ax= [rojdr+o, 


que pode ser reescrita como 


/ s(y)y' dx = / f(x)dx +. 


Fazendo a substituição y'dx = dy, obtemos 


femáy= f fixe. 


Atenção: Não se deve memorizar a fórmula obtida no final. O que fize- 
mos aqui foi mostrar o caminho que deve ser seguido para resolver uma 
equação separável. 


As curvas, que são soluções de uma equação separáve,l podem ser vistas 
como curvas de nível da função 


z = F(x,y) =h(y(5))) — f fx. 


Exemplo 2.1. Vamos, agora, encontrar a solução geral da equação diferen- 
cial 


e =—4x ou 2yy = —4x. 


Integrando-se em relação a x ambos os membros, obtemos 


foy dx = - farir+c. 


Fazendo a substituição y'dx = dy, obtemos 
[way = — faxdr+c. 


Assim, a solução geral é dada implicitamente por 
= -2x +C. 
As soluções são elipses (Figura 2.1), que são as curvas de nível da função 
z = F(x,y) = 2 +22. 


O gráfico da função F(x,y) = y? + 2x? é um paraboloide elíptico (Figura 
2.2). 


AULA 2 


Figura 2.1 - Soluções da equação diferencial do Exemplo 2.1 


Figura 2.2 - Soluções da equação diferencial do Exemplo 2.1 
como curvas de nível do paraboloide elíptico z = F(x,y) = 2x? + y7 


Exemplo 2.2. (a) Encontre a solução do problema de valor inicial 


dy _ 2x-1 
| dx 3-3 


y(1) = 0. 


(b) Determine o intervalo de validade da solução, ou seja, o maior inter- 
d 
valo contendo xo = 1 para o qual a solução y(x) e sua derivada = 


estão definidas. 
(c) Determine os pontos onde a solução tem um máximo local. 


(d) Faça um esboço do gráfico da solução. 
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Solução: 


(a) Podemos reescrever a equação como 


(c 


(d 


<~“ 


x 


~ 


(3y? —3)y' = 2x — 1. 
Integrando-se em relação a x ambos os membros, obtemos 


Jer — 3)y' dx = fo = Där+ C. 
Fazendo a substituição y'dx = dy, obtemos 

fer -“3)dy= f(x -1)dx+C. 
Assim, a solução geral é dada implicitamente por 


y = 3y =x? -x +C. 


Para encontrar a solução que satisfaz a condição inicial y(1) = 0, 
substituímos x = 1e y = 0 na solução geral, obtendo C = 0. Assim, 
a solução do problema de valor inicial é dada implicitamente por 


yY -3y -x +x =0. 


Para determinar o intervalo de validade da solução do PVI, vamos 
determinar o maior intervalo que contém x = 1, em que a solução e 


i a usa a. dy 2x-1 
sua derivada estão definidas. Pela equação + = 3223” temos que os 


pontos onde a derivada não está definida são aqueles tais que 3y? — 


3 = 0, ou seja, y = +1. Como o ponto inicial é (1,0), então a solução 
do PVI está contida na região do plano —1 < y < 1. Substituindo-se 
y = —1 na equação que define a solução, obtemos a equação x? — 
x—2 = 0, que tem solução x = —1 e x = 2. Substituindo-se y = 1 


na equação que define a solução y? — 3y — x? + x = 0, obtemos a 
equação x? — x +2 = 0, que não tem solução real. 


Como a solução está definida para todo x, mas a derivada não está 
definida para x = —1 e x = 2 e o ponto inicial xy = 1 está entre 
os valores x = —1 e x = 2, concluímos que o intervalo de validade 
da solução é o intervalo (—1,2), que é o maior intervalo em que a 
solução y(x) e a sua derivada estão definidas. 


Nos pontos onde a solução tem máximo local, a reta tangente à curva 
é horizontal, ou seja, pontos onde Ed = 0. Neste caso não precisa- 
mos calcular a derivada da solução, pois a derivada já está dada pela 
equação diferencial, ou seja, 

dy  2x—1 


dx 3-3 


Assim, a reta tangente é horizontal para x tal que 2x — 1 = 0, ou seja, 
somente para x = 1/2. 


Nos pontos x = —1 e x = 2, a reta tangente à curva solução y? — 3y — 
x? + x = 0 é vertical, ou seja, E = 0, pois, pela equação diferencial, 


dx 1 3-3 


dy E 2x — 1" 


AULA 2 


para x £ 1/2. Assim, já sabemos, pelo item (b), que a solução está 
contida em uma curva que passa pelos pontos (—1,—1) e (2,—1), 
onde a tangente é vertical, e que passa pelo ponto inicial (1,0). Neste 
ponto a inclinação da tangente é —1/3, pois, substituindo-se x = 1 


ey = 0 na equação diferencial, obtemos Ed = —1/3. Além disso, 
sabemos que o único ponto em que a tangente é horizontal ocorre 
para x = 1/2. Deduzimos daí que a solução é crescente até x = 1/2, 
depois começa a decrescer. 


Figura 2.4 - Soluções da equação diferencial e do problema de valor inicial 
do Exemplo 2.2 
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Figura 2.5 - Soluções da equação diferencial do Exemplo 2.2 como curvas 
de nível de uma função de duas variáveis z = f (x,y) = y? — 3y — x? + x 
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AULA 2 


1. Resolva as equações: 
fa o O 
(b) y°- 1- (2y + xy)y' =0. 
(c) (ayx? + by)y' — x = 0 para a,b € R,a £0. 
(d) (ax? + b)!/?y' — xy? = 0 para a,b € Ra £ 0. 
(e) (ay? +b)! — xyy' = 0 para a,b € Ra £ 0. 
(£) ay? +b — x?yy' = 0 para a,b € R,a £ 0. 


2. (a) Encontre a solução do problema de valor inicial 


dy _ 2x+1 
dx 342-3 


y(0) = 0. 


(b) Determine o intervalo de validade da solução. 
(c) Determine os pontos onde a solução tem um máximo local. 
(d) Faça um esboço do gráfico da solução. 


3. Mostre que a equação linear y' + p(t)y = q(t) é equivalente a uma equação separável se 
(a) p(t)=aeg(t)=b,paraa,bER. 
(b) p(t) = q(t). 
(c) q(t) =0. 
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2.2 Aplicações 


OBJETIVOS 
Ao terminar esta seção você deverá ser capaz de: 


e Modelar várias aplicações que levam a equações diferenciais de 1º ordem. 
e Identificar o tipo da equação diferencial correspondente ao modelo obtido. 


e Resolver um problema de valor inicial correspondente a uma equação diferencial 
obtida no modelo. 


2.2.1 Dinâmica Populacional 
Crescimento Exponencial 


O modelo mais simples de crescimento populacional é aquele em que se 
supõe que a taxa de crescimento de uma população dy é proporcional à 
população presente naquele instante y(t). Podemos descrever o problema 
de encontrar y(t) como o problema de valor inicial 


dy 
a“ 
y(0) = yo. 


A equação é linear e pode ser reescrita como 
= —ky=0. (2.4) 
Para resolvê-la vamos determinar o fator integrante 
— pf—kdt _ „—kt 
u(t) = e = e". 


Multiplicando-se a equação (2.4) por u(t) = e*t, obtemos 


gey) =0. 
Integrando-se ambos os membros, obtemos 
eyt) =C ou y(t) = Ch. 
Substituindo-se t = 0 e y = yọ, obtemos 
yo = CA? = C. 


Ou seja, a solução do problema de valor inicial é 


y(t) = yoe". 


Exemplo 2.3. Consideremos uma situação formada por uma população 
de organismos zooplanctônicos. São colocadas em um béquer 3 fêmeas 
partenogenéticas grávidas (não há necessidade de fecundação pelo ma- 
cho) de um microcrustáceo chamado cladócero em condições ideais de 
alimentação, temperatura, aeração e iluminação e ausência de predadores. 
Sabendo-se que em 10 dias havia 240 indivíduos, determine a população 
em função do tempo supondo-se que a taxa de crescimento da população 
é proporcional à população atual (crescimento exponencial). 


AULA 2 


A população, y(t), é a solução do problema de valor inicial 


dy 
ri 
y(0)=3 


que, como vimos acima, tem solução 
kt kt 
y(t) = yoe” = 3e“. 


Como em 10 dias a população é de 240 indivíduos, então, substituindo-se 
t = 10 e y = 240, obtemos 


In 80 
240 = 3e!0 k= 
0 = 3e > 10 
Assim, a função que descreve como a população de bactérias varia com o 
tempo é 


oop esce ENSENEREESECENEENEEETEES o opd 


100 F 


—100 
-5 


Figura 2.6 - Solução do problema do Exemplo 2.3 e dados obtidos experi- 
mentalmente 
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Tabela 2.1 - Número de indivíduos por litro de uma população de cladóceros 
(Daphnia laevis)em experimento de laboratório (dados obtidos de [3]) 


Dias | População || Dias | População 
1 3 13 510 
2 7 14 630 
3 10 15 638 
4 9 16 628 
5 39 17 666 
6 39 18 668 
7 40 19 620 
8 113 20 663 
9 180 21 667 
10 240 22 645 
11 390 23 690 
12 480 24 650 


Crescimento Logístico 


Para levar em conta que a população y(t) tem um valor máximo sustentável 
ym, podemos supor que a taxa de crescimento, além de ser proporcio- 
nal à população atual, é proporcional também à diferença entre ym e à 
população presente. Neste caso a população como função do tempo, y(t), 
é a solução do problema de valor inicial 


A equação é separável. Multiplicando-se a equação por 7 a a) obtemos, 


1 
yum- y)” 


Integrando-se em relação a t, obtemos 


1 , 
dt= f kdt+ Ci. 
Err 1 


Fazendo-se a substituição y'dt = dy, obtemos 


1 
aiy = [kát+ Cr 
| au” i 


Para calcular a integral do lado esquerdo vamos decompor 


t=k (2.5) 


ylym—=y) 
frações parciais 
1 A B 


yym-y) Y ym—y 
Multiplicando-se a equação acima por y(ym — y), obtemos 
1 = A(ym — y) + By. 
Substituindo-se y = 0 e y = ym, obtemos A = 1/yme B = 1/ym. Assim, 


1 1/71 1 1 
ii Í dy+ f dy) = nu-i . 
J au n” YM ( 1) ay!) a e 


AULA 2 


Logo, a solução da equação é dada implicitamente por 
In|y| = In |ym -y| = kymt + Gi. 
Usando propriedades do logaritmo, podemos reescrever como 


y 
YMY 


In 


| = C1 F kymt. 


Aplicando a exponencial a ambos os membros e eliminando-se o valor 
absoluto, obtemos 


Uvas +eCreymkt — Coymkt, 
yYM—Y 


Observe que, como C; é uma constante, então +el1 também é uma cons- 
tante, que chamamos de C. Substituindo-se t = tọ e y = yo na equação 


acima, obtemos 
y 0 ey M k to 


C e 
YM — Yo 
Vamos explicitar y(t). 
v=(ym-y)Ceomt > y+Comtty = yyCetmtt, 
Portanto, a solução do problema de valor inicial é 


kt Yoym oymk(t—to) k(t-t 
= Cyme!M = Mo? yoy me™ ( 0) 


1+Ceymkt 14 epo) Ym — Yo + yocymklt-to)" 


y(t) 


Dividindo-se numerador e denominador por eYmkt, obtemos 


YoyYM 
vo + (ym — yo)eYmklt-to) | 


y(t) = 


Observe que 
lim y(t) = ym. 


t— oo 


Exemplo 2.4. Consideremos a mesma situação do Exemplo 2.3, ou seja, 
são colocadas em um béquer 3 fêmeas partenogenéticas grávidas (não há 
necessidade de fecundação pelo macho) de um microcrustáceo chamado 
cladócero em condições ideais de alimentação, temperatura, aeração e 
iluminação e ausência de predadores. Sabendo-se que essa população 
atinge o máximo de 690 indivíduos e que em 10 dias havia 240 indivíduos, 
determine a população em função do tempo, supondo-se que a taxa de 
crescimento da população é proporcional tanto à população atual quanto 
à diferença entre a população máxima e a população atual (crescimento 
logístico). 

A população como função do tempo, y(t), é a solução do problema de valor 
inicial 


Yo ky (690 — y) 


dt 
y(0) = 3. 
A equação é separável. Multiplicando-se a equação por T’ obtemos 
: =k (2.6) 


y(690 =)” T+ 
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Integrando-se em relação a t, obtemos 


1 I 
dt= [kdt+C. 
| y(690 — y)” / rs 


Fazendo-se a substituição y'dt = dy, obtemos 


1 
= y= f ht 
len Y / E 


Para calcular a integral do lado esquerdo, vamos decompor TEN em 


frações parciais 
1 A B 


y(690 —y) y 690-y 
Multiplicando-se a equação acima por y(690 — y), obtemos 


1 = A (690 — y) + By. 
Substituindo-se y = 0 e y = 690, obtemos A = 1/690 e B = 1/690. Assim, 


1 1 1 1 1 
| zomp” = (J 3w | [ mst) szg (in lvl = In [690 — y|). 


Logo, a equação (2.6) tem solução dada implicitamente por 


In |y| — In |690 — y| = k690t + C1. 


Usando propriedades do logaritmo, podemos reescrever como 


y 
In |—=— | = Cy + k690t. 
E | 690 — | ai 
Aplicando a exponencial a ambos os membros, obtemos 
Y peço — Ce690kt 
690 +e~te Re 


Observe que, como Cj é uma constante, então +eQ1 também é uma cons- 
tante, que chamamos de C. Substituindo-se t = 0 e y = 3 na equação 
acima, obtemos 

3 3 1 


C= 0—3 T 6877 229 


Vamos explicitar y(t). 
y = (690 — y)Ce6™ => y+ Cey = 690C e, 


Portanto, a solução do problema de valor inicial é 


690C e690kt 690e690kt 690e690kt 
y(t) = 14 Ce690kt — 1/C + e690kt — 229 + e690kt” 
—690kt 


ou, multiplicando-se numerador e denominador por e , obtemos como 
solução do PVI 
moe 690 l 
229e-690kt + 1 
Para determinar o valor de k, vamos usar o fato de que em 10 dias havia 
240 indivíduos. Substituindo-se t = 10 e y = 240, obtemos 
690 690 23 15 


— 22 —6900k _ 1= 1= 
229e-6900F + 1 do 240 8 8 


(2.7) 


240 


AULA 2 


In a 
> —690k = a 


Logo, substituindo-se o valor de —690k obtido acima na solução do PVI 
(2.7), obtemos que a população de cladóceros em função do tempo é dada 


por 
690 690 


y(t) = 15 = T 


In 


229º 1! +1 229 (rãs) 


10 


+1 


Figura 2.7 - Solução do problema de valor inicial do Exemplo 2.4 e dados 
obtidos experimentalmente 
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2.2.2 Datação por Carbono 14 


A proporção de carbono 14 (radioativo), em relação ao carbono 12 presente 
nos seres vivos, é constante. Quando um organismo morre a absorção de 
carbono 14 cessa, e a partir de então o carbono 14 vai se transformando em 
carbono 12 a uma taxa que é proporcional à quantidade presente. Podemos 
descrever o problema de encontrar a quantidade de carbono 14 em função 
do tempo, y(t), como o problema de valor inicial 


dy _ 
e —ky. 
y(0) = yo. 


A equação é a mesma do crescimento exponencial (trocando-se k por —k), 
e vimos na página 30 que este problema tem solução 


y(t) = yoe ™, 


em que yọ é a quantidade no instante t = 0. 


Exemplo 2.5. Em um pedaço de madeira é encontrado 1/500 da quan- 
tidade original de carbono 14. Sabe-se que a meia-vida do carbono 14 é 
de 5600 anos, ou seja, que em 5600 anos metade do carbono 14 presente 
transformou-se em carbono 12. Vamos determinar a idade deste pedaço 
de madeira. 


O problema de valor inicial que descreve esta situação é 


dy 
LAN, 
dt Y 
y(0) = yo, 
que tem solução 
y(t) = yoe ™. 
Substituindo-se t = 5600 e y = yo/2 (meia-vida), obtemos 
In2 
— 4,0 k:5600 2 
Yo/2 = yoe > k= 00 
Agora, substituindo-se y = yo /500, obtemos 
vo -kt In500 5600 1n 500 
40 — = = = 502 ; 
500 7 “0º t k n? 50200 anos 
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Yo 


Yo!2 pinapan N 


5000 10000 15000 20000 25000 30000 35000 40000 45000 50000 


Figura 2.8 - Solução do problema de valor inicial do Exemplo 2.5 


2.2.3 Misturas 


Figura 2.9 - Tanque 


Vamos supor que um tanque contenha uma mistura de água e sal com um 
volume inicial de Vo litros e Qọ gramas de sal, e que uma solução salina 
seja bombeada para dentro do tanque a uma taxa de T litros por minuto 
possuindo uma concentração de Ce gramas de sal por litro. Suponha que a 
solução bem misturada saia a uma taxa de T; litros por minuto. 


A taxa de variação da quantidade de sal no tanque é igual à taxa com que 
entra sal no tanque menos a taxa com que sai sal do tanque. 


A taxa com que entra sal no tanque é igual à taxa com que entra a mis- 
tura, Te, vezes a concentração de entrada, Ce. E a taxa com que sai sal 
do tanque é igual à taxa com que sai a mistura do tanque, Ts, vezes a 
concentração de sal que sai do tanque, Cs. Como a solução é bem mis- 
turada esta concentração é igual a concentração de sal no tanque, ou seja, 


Colt) = oi) E 
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Como o volume no tanque, V(t), é igual ao volume inicial, Vo, somado ao 
volume que entra no tanque menos o volume que sai do tanque, então 


V(t) = Vo + Tet — Tot = Vo + (Te — To)t. 


Assim, a quantidade de sal no tanque, Q(t), é a solução do problema de 
valor inicial 


do Q 
e RCE “Va + (T; Tot 
Q(0) = Qo 


Exemplo 2.6. Num tanque há 100 litros de salmoura contendo 30 gramas 
de sal em solução. Água (sem sal) entra no tanque à razão de 6 litros por 
minuto e a mistura se escoa à razão de 4 litros por minuto, conservando-se 
a concentração uniforme por agitação. Vamos determinar qual a concen- 
tração de sal no tanque ao fim de 50 minutos. 


O problema pode ser modelado pelo seguinte problema de valor inicial 
io 4 Q 
dt 100 + 2t 
Q(0) = 30. 
A equação é linear e pode ser escrita como 
dQ Q 
H4 = 
dt 100 + 2t j 
Um fator integrante é, neste caso, 


u(t) = ef morat — e21n(100+2t) — pln((100+2t)?°) — (100 +24). 
Multiplicando-se a equação por u(t) = ef wrat — (100 + 2t)?, obtemos 


ç ((100+220) =0. 


Integrando-se, obtemos 

(100 +24 Q(t) =C, 
ou seja, é 

Q(t) = (100 4207" 
Substituindo t = 0 e Q = 30, obtemos 

C = 30 - 100° = 3 - 10°. 
Substituindo o valor de C encontrado, obtemos 
3-105 

(100 + 2t)?" 


A concentração é o quociente da quantidade de sal pelo volume, que é 
iguala V(t) = 100 + 2t. Assim, 


Q(t) = 


3-105 
tj = = 
e(O = 00 +23 
e após 50 minutos 
3:10 3 


“(50) = Ta00j3 — 80 


= 0,0375 gramas /litro. 
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Figura 2.10 - Solução do problema de valor inicial do Exemplo 2.6 


Figura 2.11 - Concentração como função do tempo para o problema do 
Exemplo 2.6 
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2.2.4 Lei de Resfriamento de Newton 


A lei de resfriamento de Newton diz que a taxa de variação da tempera- 
tura T(t) de um corpo em resfriamento é proporcional à diferença entre a 
temperatura atual do corpo T(t) e a temperatura constante do meio ambi- 
ente Tm, ou seja, a temperatura do corpo, T(t) é a solução do problema de 
valor inicial 


dT 
| AE = k(T == Tm) 


T(0) = To. 


Exemplo 2.7. O café está a 90º C logo depois de coado e, um minuto de- 
pois, passa para 85º C, em uma cozinha a 25º C. Vamos determinar a tem- 
peratura do café em função do tempo e o tempo que levará para o café 
chegar a 60º C. 


T(0) = 90, T(1) = 85. 
Dividindo-se a equação por T — 25 


dT 
| 2T =HT-25) 


Integrando-se em relação a t 


1 
T'dt = Í kdt, 
J T=25 


[rr = [ras 
In|T-25]=kH+G, 
T(t) = 25 + ele = 25 + Cet 
Substituindo t = 0 e T = 90 


9 =25+C > C=65. 
T(t) = 25 + 65e, 
Substituindo-se t = 1 e T = 85 
85-25 +65 = k=In(60/65). 
Assim, a temperatura do café em função do tempo é dada por 
T(t) = 25 + 65¢™(60/65)t, 


Substituindo T = 60 


Logo, o tempo necessário para que o café atinja 60° é de 


= n(35/65) 
= In(60/65) 


x 8min. 
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Figura 2.12 - Solução do problema de valor inicial do Exemplo 2.7 


2.2.5 Juros 


Vamos supor que façamos uma aplicação de uma quantia Sọ em um banco 
e que a taxa de variação do investimento dS seja proporcional ao saldo em 
cada instante S(t). Podemos descrever o problema de encontrar S(t) como 


o problema de valor inicial 
dS 
E rS. 
50) = So. 


Este problema já resolvemos antes e tem solução 
S(t) = Soe”. (2.8) 


Pode parecer que este modelo não seja muito realista, pois normalmente 
os juros são creditados em períodos inteiros igualmente espaçados. Ou 
seja, se j é a taxa de juros em uma unidade de tempo, então o saldo após n 
unidades de tempo S(n) é dado por 


S(1) = So+Soj = So(1 +j) 
s(2) = S(1)(1+j) = So(1 + j}? 


S(n—1)(1 +j) = So(1 + j)”. (2.9) 


Substituindo-se t por n na solução do problema de valor inicial obtida (2.8) 
e comparando com (2.9), obtemos que 


Soe” = So(1 Po 


S(n) 


ou seja, 
1+j=" ou r=ln(1+j). (2.10) 


Assim, a hipótese inicial de que os juros são creditados continuamente é 
realista desde que a constante de proporcionalidade na equação diferen- 
cial r e a taxa de juros j estejam relacionadas por (2.10). Para pequenas 
taxas de juros os dois valores são muito próximos. Por exemplo, j = 4% 
corresponde a r = 3,9 % e j = 1% corresponde a r = 1%. 
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Exemplo 2.8. Vamos supor que uma aplicação renda juros de 1% ao mês 
(continuamente). Vamos encontrar o saldo como função do tempo e o saldo 
após 12 meses se o saldo inicial é de R$ 100, 00. 


Podemos descrever o problema de encontrar S(t) como o problema de va- 
lor inicial 


ds 

> =0,01 
q —001S 
S(0) = 100. 


Este problema já resolvemos antes e tem solução 
S(t) = 100e, 


Assim, em 12 meses o saldo é 


S(12) = 100€20112 ~ R$ 112,75. 


114 


112} 


1410F 


Saldo em R$ 


Meses 


Figura 2.14 - Saldo em função do tempo para o problema do Exemplo 2.8 
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Vamos supor, agora, que, além do investimento inicial So, façamos depó- 
sitos ou saques continuamente a uma taxa constante d (positivo no caso de 
depósitos e negativo no caso de saques), então, neste caso, o modelo que 
descreve esta situação é o do problema de valor inicial 


dS 
E d 
Er rS + 
S(0) = So. 
A equação é linear e pode ser reescrita como 
dS 
— -rS =d. 2.11 
m S (2.11) 


Para resolvê-la vamos determinar o fator integrante 
u(t) Z eJ —rdt — et, 
Multiplicando-se a equação (2.11) por u(t) = e7", obtemos 


d 
ate "S) = de". 


Integrando-se ambos os membros, obtemos 
d d 
e™™ S(t) = Za +C ou S(t) = Ce" — E 


Substituindo-se t = 0 e S = So, obtemos 


o d d 
So = Ce"? + Debi 
r r 
Ou seja, a solução do problema de valor inicial é 
d 
S(t) = Soe + z(e” —1). (2.12) 


Vamos comparar este resultado com o caso em que, além dos juros serem 
creditados em intervalos constantes, os depósitos ou saques de valor D 
são feitos em intervalos constantes. Neste caso o saldo após n unidades de 
tempo é dado por 


S) = So(1+j)+D 

S(2) = SU+HIP+DA+N)+D 

S(n) = So1+j"+D((1 Fj E e i) 

Sa) = WD = é: (2.13) 


Foi usada a soma de uma progressão geométrica. Substituindo-se t por 
n na solução do problema de valor inicial (2.12) e comparando-se com a 
equação (2.13), obtemos que 
d 1+/)"-1 
Soe™ + (e —1) = So(1 | j)” | pt + a 


ou seja, 
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Usando (2.10), obtemos 
(e"— 1)d 


In(1+7)D 
E DO a pe E (2.14) 
j r 
Assim, podemos também neste caso usar o modelo contínuo em que os 
depósitos ou saques são feitos continuamente desde que a taxa contínua 
de depósitos d e os depósitos constantes D estejam relacionados por (2.14). 


d 


Exemplo 2.9. Suponha que seja aberta uma caderneta de poupança com o 
objetivo de no futuro adquirir um bem no valor de R$ 40.000, 00. Suponha 
que os juros sejam creditados continuamente a uma taxa de r = 1% ao mês 
e que os depósitos também sejam feitos continuamente a uma taxa cons- 
tante, sendo no início o saldo igual a zero. Vamos determinar de quanto 
deve ser a taxa de depósito mensal para que em 20 meses consiga atingir o 
valor pretendido. 


ds 1 

dt 100 

S(0) — 0. 
A equação é linear e pode ser reescrita como 


ds 1 
r mo = É (2.15) 


Para resolvê-la precisamos determinar o fator integrante 


S+d 


u(t) = ef -4t = e- mt 


Multiplicando-se a equação (2.15) por u(t) = e mf, obtemos 


E 
dt 


Integrando-se ambos os membros, obtemos 


(e 100tS) = de mt, 


em! S(t) = —100de 10! +C ou S(t) = Certo! — 1004. 
Substituindo-se t = 0 e S = 0, obtemos 
0=Cermº-100d => C=100d. 
Ou seja, a solução do problema de valor inicial é 
S(t) = 100d(e1o! — 1). (2.16) 
Substituindo-se t = 20 e S = 40000, obtemos 


40000 = 1004 (eñ — 1), 


400 400 


d = N 
eo —1 0,22 


x R$ 1818,18. 


Esta é a taxa de depósito mensal, supondo-se que os depósitos sejam rea- 
lizados continuamente. Vamos determinar o depósito mensal correspon- 
dente. zdi 
r—1i)jd ^1 — 1)1818,1 
pe ta 81818 pg1827,30. 
r 0,01 
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Figura 2.15 - Saldo em função do tempo quando são feitos depósitos a uma 


taxa constante 


Figura 2.16 - Solução do problema de valor inicial do Exemplo 2.9 


= 
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. Um tanque contém inicialmente 100 litros de água pura. Então, água salgada, contendo 30 e” 1 


Exercícios 


. Um tanque contém 100 litros de uma solução a uma concentração de 1 grama por litro. Uma 


= a =. : 
solução com uma concentração de 2te” 10! gramas por litro entra no tanque a uma taxa cons- 
tante de 1 litro por minuto, enquanto a solução bem misturada sai à mesma taxa. 


(a) Determine a quantidade de sal no tanque em cada instante t, onde t é contado a partir do 
início do processo. 


(b) Calcule a concentração de sal no tanque t = 10 minutos após o início do processo. 


t 


gramas de sal por litro, passa a ser bombeada para o tanque a uma taxa de 10 litros por minuto. 
Simultaneamente a solução passa a ser agitada e retirada do tanque na mesma taxa. 


(a) Determine a quantidade de sal no tanque em cada instante t, onde t é contado a partir do 
início do processo. 


(b) Calcule em que instante a concentração de sal no tanque será de 7,5 gramas por litro. 


. Um tanque contém inicialmente 100 litros de água e 100 gramas de sal. Então, uma mistura de 


água e sal na concentração de 5 gramas de sal por litro é bombeada para o tanque a uma taxa 
de 4 litros por minuto. Simultaneamente a solução (bem misturada) é retirada do tanque na 
mesma taxa. 


(a) Determine a quantidade de sal no tanque em cada instante t, onde t é contado a partir do 
início do processo. 


b) Calcule a concentração limite de sal no tanque quando t — œ e o tempo necessário para 
que q P P 
que a concentração atinja metade deste valor. 


. Suponha que um tanque contenha uma mistura de água e sal com um volume inicial 100 litros 


e 10 gramas de sal e que uma solução salina seja bombeada para dentro do tanque a uma taxa 
de 3 litros por minuto possuindo uma concentração de 1 grama de sal por litro. Suponha que a 
solução bem misturada saia a uma taxa de 2 litros por minuto. 


(a) Determine a quantidade de sal no tanque em cada instante t, onde t é contado a partir do 
início do processo. 


(b) De qual valor se aproxima a concentração quando o tanque está enchendo, se a sua capaci- 
dade é de 200 litros? 


. Suponha que um tanque contenha uma mistura de água e sal com um volume inicial 100 litros 


e 10 gramas de sal e que água pura seja bombeada para dentro do tanque a uma taxa de 1 litro 
por minuto. Suponha que a solução bem misturada saia a uma taxa de 2 litros por minuto. 


(a) Determine a quantidade de sal no tanque em cada instante t, onde t é contado a partir do 
início do processo. 


(b) De qual valor se aproxima a concentração quando o tanque se aproxima de ficar vazio? 


. Num processo químico, uma substância se transforma em outra, a uma taxa proporcional à 


quantidade de substância não transformada. Se esta quantidade é 48 ao fim de 1 hora, e 27, ao 
fim de 3 horas, qual a quantidade inicial da substância? 


. A população de bactérias em uma cultura cresce a uma taxa proporcional ao número de bactérias 


no instante t. Após três horas, observou-se a existência de 400 bactérias. Após 9 horas, 2500 
bactérias. Qual era o número inicial de bactérias? 
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Uma população de bactérias cresce a uma taxa proporcional à população presente. Sabendo-se 
que após uma hora a população é 2 vezes a população inicial, determine a população como 
função do tempo e o tempo necessário para que a população triplique. Faça um esboço do 
gráfico da população em função do tempo. 


Suponha que em uma comunidade de 100 pessoas inicialmente apenas uma pessoa seja por- 
tador de um vírus e que a taxa com que o vírus se espalha na comunidade seja proporcional 
tanto ao número de pessoas infectadas como também ao número de pessoas não infectadas. Se 
for observado que após 4 semanas 5 pessoas estão infectadas. Determine o número de pessoas 
infectadas em função do tempo. Faça um esboço do gráfico da solução. 


Um termômetro é levado, de uma sala onde a temperatura é de 20° C, para fora onde a tempe- 
ratura é de 5° C. Após 1/2 minuto o termômetro marca 15° C. 

(a) Determine a temperatura marcada no termômetro como função do tempo. 

(b) Qual será a leitura do termômetro após 1 minuto? 

(c) Em quanto tempo o termômetro irá marcar 10° C? 


Suponha que um automóvel sofra depreciação continuamente numa taxa que é proporcional ao 
seu valor num instante t. Este automóvel novo custa R$ 35000,00. Após um ano de uso o seu 
valor é R$ 30000,00. Qual será o valor do automóvel após dois anos de uso? 


Com o objetivo de fazer uma previdência particular uma pessoa deposita uma quantia de R$ 

100,00 por mês durante 20 anos (suponha que o depósito seja feito continuamente a uma taxa 

de R$ 100,00 por mês). 

(a) Supondo que neste período a taxa de juros seja de 1 % ao mês (contínua), qual o valor que 
esta pessoa iria ter ao fim deste período. 

(b) Se após o período anterior esta pessoa quisesse fazer retiradas mensais, qual deveria ser o 
valor destas retiradas para que em 20 anos tenha desaparecido o capital, se a taxa de juros 
continuasse em 1 % (contínua)? 
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Equações Diferenciais de 1º Ordem 
3º Parte 


3.1 Análise Qualitativa 


3.1.1 Equações Autônomas 


As equações autônomas são equações da forma 
dy 
“Ed (y). (3.1) 


Vamos supor que f (y) seja derivável com derivada contínua no intervalo 
de estudo. Para as equações autônomas podemos esboçar várias soluções 
sem ter que resolver a equação, pois a equação diferencial fornece a incli- 


nação da reta tangente às soluções, z, como função de y e assim pode- 

mos saber como varia com y o crescimento e o decrescimento das soluções. 

Além disso, podemos saber os valores de y para os quais as soluções têm 

pontos de inflexão e como varia a concavidade das soluções com y, pois 
dy ddy d 


a2 = iat E a O) 


e pela Regra da Cadeia 


EF = FN = Ff). 


Assim, g 
=. 


Observe que se y1,...,Yyy são zeros da função f(y), então y(t) = y; são 
soluções constantes da equação (3.1), para i = 1,...,k (verifique!). 


Definição 3.1. (a) Sejam y1,..., Yg zeros da função f(y). Os pontos y; 
são chamados pontos críticos ou de equilíbrio da equação (3.1), e 
as soluções y(t) = y; são denominados soluções de equilíbrio ou 
estacionárias da equação (3.1). 
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(b) Um ponto de equilíbrio y; é chamado estável se para y(to) um pouco 
diferente de y;, y(t) se aproxima de y;, quando t cresce. 


(c) Um ponto de equilíbrio y; é chamado instável se para y(to) um pouco 
diferente de y;, y(t) se afasta de y;, quando t cresce. 


d 
Figura 3.1 - = = f(y) nas proximidades de um ponto de equilíbrio estável 


d 
Figura 3.2 - Soluções de Z = f(y) nas proximidades de um ponto de 


equilíbrio estável 


O ponto de equilíbrio y; é estável se f(y) < 0 para y próximo de y; com 
y > yie f(y) > 0 para para y próximo de y; com y < yi. Pois neste caso 


e se y(to) é um pouco maior do que y;, então a derivada Y = f(y) 
é negativa e, portanto, a solução y(t) é decrescente e assim y(t) se 
aproxima de y;, quando t cresce. 


e se y(to) é um pouco menor do que y;, então a derivada Y = f(y) é 
positiva e, portanto, a solução y(t) é crescente e assim y(t) se apro- 
xima de y;, quando t cresce. 


O ponto de equilíbrio y; é instável se f(y) > 0 para y próximo de y; com 
y > yie f(y) < 0 para para y próximo de y; com y < yi. Pois neste caso 


e se y(to) é um pouco maior do que y;, então a derivada A = f(y)é 


positiva e, portanto, a solução y(t) é crescente e assim y(t) se afasta 
de y;, quando t cresce. 
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d 
Figura 3.3 - a = f(y) nas proximidades de um ponto de equilíbrio 
instável 
v=y(0) 
yi | 
ESSE 
: . dy o 
Figura 3.4 - Soluções de Ta f(y) nas proximidades de um ponto de 


equilíbrio instável 


d 
e se y(to) é um pouco menor do que y;, então a derivada = = f(y) 
é negativa e, portanto, a solução y(t) é decrescente e assim y(t) se 
afasta de y;, quando t cresce. 


Exemplo 3.1. Considere a equação diferencial: 


dy do 
g TY (3.2) 


Vamos esboçar várias soluções da equação. Para isto vamos obter os pon- 
tos de equilíbrio. Depois vamos determinar como varia o crescimento e o 


decrescimento das soluções com y. E finalmente para quais valores de y as 
soluções têm ponto de inflexão. 


Os pontos de equilíbrio são as raízes de y? — y = 0, ou seja, y1 = 0ey = 1. 


d 
Como c =y? —y < 0, para 0 < y < 1, então as soluções são decrescentes 


para0<y<1. 


d 
Como A =7 -y>0, para y < 0 e para y > 1, então as soluções são 
crescentes para y < 0 e para y > 1. 
Vamos determinar para quais valores de y as soluções têm pontos de in- 
flexão e como varia a concavidade das soluções com y calculando a se- 
gunda derivada. 
dy ddy d ( 
de Catar AY 


y). 
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Mas, pela Regra da Cadeia, 


de? = dy _ 2 
q T») = Quo Ds = (Quo yo y). 
Assim, 
d? 
=a = (Cy D(y -y). 
Logo, as soluções têm pontos de inflexão para y = 1/2, y = 0e y = 1. 


Observamos que o ponto de equilíbrio y = O é estável, pois para valores 
de y próximos de y1 = 0 as soluções correspondentes y(t) estão se aproxi- 
mando de y = 0, quando t cresce. O ponto de equilíbrio y2 = 1 é instável, 
pois para valores de y próximos de y2 = 1 as soluções correspondentes 
y(t) estão se afastando de y2 = 1, quando t cresce. Com as informações 
sobre os pontos críticos, regiões de crescimento e decrescimento, pontos de 
inflexão, podemos fazer um esboço dos gráficos de algumas soluções. 


Figura 3.5 - a = f(y) da equação 3.2 
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Figura 3.6 - Algumas soluções da equação 3.2 


3.1.2 Campo de Direções 


Uma maneira de se ter uma ideia do comportamento das soluções de uma 
equação diferencial de 1º ordem 


dy 
I — f(t 
dt f( 1 y), 
sem ter de resolvê-la, é desenhar o campo de direções 
1 dy 1 


ED Ra) = Tt) 


da seguinte forma: 


(a) Constrói-se uma malha retangular consistindo em pelo menos uma 
centena de pontos igualmente espaçados; 


(b) Em cada ponto da malha desenha-se um segmento orientado unitário 
que tem inclinação igual à da reta tangente à solução da equação que 
pelo ponto da malha, ou seja, na direção e sentido de 


dy 
(2) = (L fly) 
e com comprimento igual a 1. 


Desenhar o campo de direções é, como está dito em [?], “uma tarefa para 
a qual o computador é particularmente apropriado e você deve, em geral, 
usar o computador para desenhar um campo de direções.” Por isso escre- 
vemos uma função para o MATLABÊ que está no pacote GAAL e que torna 
esta tarefa mais fácil, chamada campo (f , [xmin xmax], [ymin ymax]). 


Entretanto, para as equações autônomas, como as que estudamos na seção 
anterior, é fácil desenhar o campo de direções, pois as inclinações variam 
somente com y. 


Para a equação do Exemplo 3.1 está desenhado a seguir o campo de direções. 


od 


E o E A A RÃS E É 
FIA 
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da 


22500, 


A O MD O A A R, 


ha 
NR 
iate 
NNSS 
RARE ap nd ses 
Eesti 
ae 
a SE E 
ANSERA 
AINNE 
NE rest 


Figura 3.7 - Campo de direções da equação do Exemplo 3.1 
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Exercícios 


Para as equações diferenciais autônomas dadas 


dy 
Z = fly) 
(a) Esboce o gráfico de f (y) em função de y, determine os pontos de equilíbrio e classifique cada 
um dos pontos de equilíbrio como assintoticamente estável ou instável. Justifique. 
(b) Determine como varia o crescimento das soluções com y. 
(c) Determine para quais valores de y as soluções têm pontos de inflexão. 


(d) Esboce algumas soluções da equação usando os resultados dos itens anteriores. 


dy 2 dy Z a 
L gm y. b gY y 
dy ara dy 2 
Para as equações diferenciais autônomas dadas 
W 
J fa) 


Esboce o gráfico de f(y) em função de y, determine os pontos de equilíbrio e classifique 
cada um deles como assintoticamente estável ou instável. Justifique. 


5 Y = (P-9W+9) A 
y 2 
i 7. a TOS +3y + 2). 
6. Z = (E —1)(y +4). 
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3.2 Existência e Unicidade de Soluções 


Considere novamente o problema de valor inicial 


dy _ 
y (to) = yo- 


Nem sempre este problema tem uma única solução, como mostra o próximo 
exemplo. 


—0,1 
a; -0.8 -05 -04 -02 0 0.2 0.4 0.6 0.8 1 


Figura 3.8 - Duas soluções do problema de valor inicial do Exemplo 3.2 


Exemplo 3.2. Considere o problema de valor inicial 


dy 

ar = VS 

y(0) = 0. 
Este problema tem duas soluções (verifique!) 


2 
y(t) = a parat > 0 


y2 (t) = 0. 


Observe que a equação diferencial acima é separável e é também autônoma, 
sendo y2(t) = 0 uma solução de equilíbrio. 


o) 
Se a função f(t,y) e a sua derivada o forem contínuas em um retângulo 


em torno de (to, y0), O que ocorreu no exemplo anterior não acontece como 
estabelecemos no próximo teorema. 


AULA 3 


Figura 3.9 - Retângulo em torno de (to, yo, ) onde o problema de valor ini- 
cial tem uma única solução 


Teorema 3.1 (Existência e Unicidade). Considere o problema de valor inicial 


dy _ 
y(to) = yo. 
df aait A 
Se f (t,y) e 7 são contínuas no retângulo 
R={(t,y) ER |a<t<p,8<y< y}, 


contendo (to, yo), então o problema (3.4) tem uma única solução em um intervalo 
contendo to. 


Exemplo 3.3. Para o problema de valor inicial do Exemplo 3.2, mas com o 
ponto inicial (to, 0), temos 
dy ” 
| E TVI 


y(to) = yo 


of a 
Ly) or 
fly) = vy w A 
Vemos que se (to, yo) é tal que yo > 0, então o problema de valor inicial 
acima tem solução única. 


Exemplo 3.4. Considere o problema de valor inicial 


dy _ 2 
y(to) = yo. 


Pelo Teorema 3.1, o problema de valor inicial acima tem uma única solução 
para todo (to, y0) € R2. Mas, por exemplo, para tọ = 0 e yọ = 1, o problema 


tem solução y(t) = E] (verifique!) e é válida somente no intervalo t < 1. 
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No exemplo anterior, apesar de o Teorema 3.1 garantir que em todo ponto 
(to, y0) € R? existe uma solução localmente (num intervalo em torno de 
to), estas soluções não se juntam de modo a formar soluções globais (que 
existam para todo t € R). Isto não ocorre para equações lineares, como 
provamos a seguir. 


Figura 3.10 - Solução do problema de valor inicial do Exemplo 3.4 para 
to = 0e Vo = 1 


Teorema 3.2 (Existência e Unicidade para Equações Lineares). Considere o 
problema de valor inicial 


D+ pley = a(8 
y(to) = yo. 


Se p(t) e g(t) são funções contínuas em um intervalo aberto I contendo to, então 
o problema de valor inicial tem uma única solução neste intervalo. 


Demonstração. A unicidade segue-se do Teorema 3.1 na página 61. Vamos 
provar a existência exibindo a solução do problema de valor inicial. Seja 
1 t t 
VO = to (rO) +10), emque pli) = eho, 
0 


Por hipótese a função y(t) está bem definida. Vamos mostrar que y(t) é 
solução do problema de valor inicial. 


HOYO = [ uls)alo)as+ 1 


Como p(t) e q(t) são contínuas, então 


E UOY) = KOA.: 


Derivando o produto, obtemos 


ue) + Ey = (a(o). 
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Mas ia = u(t)p(t), então a equação acima pode ser escrita como 


d 
uE) + upy = uO. 
Dividindo-se por u(t), obtemos a equação dada. 


Agora, como y(to) = yo, segue-se que y(t) dado é a solução do problema 
de valor inicial. | 


Exemplo 3.5. Considere o problema de valor inicial 


dy 2 
TR 


y(to) = yo 


2 
p(t) = 7e g(t) = t. p(t) é contínua para t Æ 0. Para ty = 2, por exem- 
plo, o problema de valor inicial tem uma única solução para t > 0, e para 
to = —3, o problema de valor inicial tem uma única solução para t < 0. 
Para tirarmos esta conclusão não é necessário resolver o problema de valor 
inicial, apesar dele estar resolvido no Exemplo 1.8 na página 21. 
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1. Determine os pontos (tọ, yọ) para os quais podemos garantir que o problema de valor inicial 


| 2 = f(t,9) 
y(to) = yo 


tem uma única solução. 


(a) Se f(t,y) = yy? — 4. (c) Se f(t, y) = y 


(b) Se f(t,y) = VFy. (d) Se f(t, y) = tyy = 1. 


2 


2. Determine o maior intervalo em que os problemas de valor inicial abaixo têm solução, sem 
resolvê-los: 


a | -DY +e-2y= z | -0% +e 
y(0) = yo y(—1) = yo 

o | Cany +ty =P E | (£ se (t+3)y = cost 
y(2) = yo y(2) = yo 
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Equações Diferenciais de 2° Ordem 
1º Parte 


4.1 Equações Homogêneas 


OBJETIVOS 
Ao terminar esta seção você deverá ser capaz de: 


e Identificar uma equação diferencial linear de 2? ordem homogênea. 

e Saber se duas funções são soluções fundamentais de uma equação diferencial 
linear de 2? ordem homogênea. 

e Dada uma solução de uma equação diferencial linear de 2º ordem homogênea, 
encontrar uma segunda solução de forma que com a primeira sejam soluções 
fundamentais da equação diferencial linear de 2º ordem homogênea. 

e Dadas duas soluções fundamentais, encontrar a solução geral de uma equação 
diferencial linear de 2º ordem homogênea. 

e Encontrar a solução geral de uma equação diferencial linear de 2° ordem 
homogênea com coeficientes constantes. 


e Resolver um problema de valor inicial correspondente a uma equação diferencial 
linear de 2° ordem homogênea. 


Uma equação diferencial linear de 2? ordem homogênea é uma equação 
da forma 2 F 
dy y = 


que pode ser reescrita como 


dy dy = 
qe TPO gr +a(by = 0 e) 


Para as equações lineares homogêneas é válido o Princípio da Superposição, 


segundo o qual se y (t) e y>(t) são soluções da equação (4.1), então 


y(t) = c1y1 (t) + coyo(t). (4.2) 


também o é para todas as constantes c; e c2. Uma expressão da forma (4.2) 
é chamada combinação linear de 11 (t) e y>(t). 


CÁLCULO V: SÉRIES NUMÉRICAS 


Vamos verificar que realmente y(t) dado por (4.2) é solução de (4.1). 


2 +p( Em +q(t)y = 
d? d 
= J2 (c1y1 + c22) 4 PD (c1y1 + coy2) + q(t) (c1y1 + coyo) 


Py Py, dyi 
ido Te cıp(t) J | cop(t) Ji cıg(t)yı + c24 (t)y2 


-o (Pn ppt, o (Bin dyz 
a (TA too aon) + e (T +O +0) 


=0:0+c,-0=0, 


pois y1 (t) e y>(t) são soluções de (4.1). Mostramos o seguinte teorema. 


Teorema 4.1 (Princípio da Superposição). Se y(t) e y>(t) são soluções de 
(4.1), então 
y(t) = cayı (t) + coyalt), 


para c1 e c constantes, também o é. 


4.1.1 Soluções Fundamentais 


Considere, agora, o problema de valor inicial 


d2y dy 
2 + AU: +a(t)y =0, 


vulto) =v0, y’ (to) = yo 


(4.3) 


em que yo e yh são condições iniciais dadas no problema. 


Vamos determinar condições sobre duas soluções y(t) e y>(t) para que 
existam constantes cj e cp tais que y(t) = cry (t) + c2y2 (t) seja solução do 
problema de valor inicial (4.3). 


Substituindo-se t = tg na solução 

y(t) = caya (t) + c2y2(t) 
e na derivada de y(t), 

y' (E) = c1y1 (t) + c22 (t), 


obtemos o sistema de equações lineares 


Eoo + coyolto) = yo 
ciy (to) + coyo(to) Yo 


que pode ser escrito na forma 


AX =B, 


66 


AULA 4 


Se a matriz do sistema A é invertível, então, para todo par de condições 
iniciais (yo, y9), O sistema tem uma única solução (c1,c>) (A solução é X = 
ATÍB). 
Mas uma matriz quadrada é invertível se, e somente se, o seu determinante 
é diferente de zero. 
Ou seja, se 
to) y2(to) 
det y1 (to 0, 
yi (to) volto) F 


então, para todo par de condições iniciais (yo, yọ), existe um único par de 
constantes (c1, c2) tal que y(t) = c1y1 (t) + c2y2(t) é solução do problema 
de valor inicial (4.3). 


Acabamos de provar o seguinte resultado. 
Teorema 4.2. Sejam y1(t) e y2(t) duas soluções da equação (4.1) tais que, em um 
ponto to € R, 
to) y2(to) 
det yı (fo 0. 
vi(to) volto) | 


Então, para todo par de condições iniciais (yo, Yo), O problema de valor inicial 


dy dy 
q2 T PUI iE +a(t)y =0, 


vlto)=yo yíto) = v) 


tem uma única solução da forma 


y(t) = c1y1 (t) + c2y2 (t). 


Definição 4.1. (a) O determinante 


to) t 
W[y1,y2] (t0) = det e lo 


é chamado wronskiano das funções y1 (t) e y>(t) em to. 


(b) Se duas soluções yı (t) e y>(t) de (4.1) são tais que o seu wronskiano 
é diferente de zero em um ponto ty, dizemos que elas são soluções 
fundamentais de (4.1). 


(c) Se yi(t) e y>(t) são soluções fundamentais de (4.1), então a família de 
soluções 


y(t) = c1y1 (t) + coyo(t), (4.4) 


para constantes c1 e c2 é chamada solução geral de (4.1). 


Assim, para encontrar a solução geral de uma equação diferencial linear 
homogênea de 2º ordem (4.1), precisamos encontrar duas soluções funda- 
mentais da equação (4.1), ou seja, duas soluções yı (t) e y>(t) tais que em 
um ponto ty € R 
det | Yilto) y2lto) | 10, 
D to) 


vilto) vol 
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Exemplo 4.1. Seja b um número real não nulo. Vamos mostrar que y1 (t) = 
cos bt e y>(t) = sen bt são soluções fundamentais da equação 


y" +by=0. 


Como y| (t) = —bsenbt, y! (t) = —b?cosbt, yh(t) = bcosbt e y4 (t) 
—b2 sen bt, então 


yl! +by = —b?° cos bt + b? cos bt = 0 


y3 + byn = —b? sen bt + b? sen bt = 0. 
Assim, y1 (t) e y>(t) são soluções da equação y” + b2y = 0. Além disso, 


| = det | no o. = b(cos? bt + sen? bt) 


=b+0 paratodot €R. 


Portanto, y(t) = cosbt e y2(t) = senbt são soluções fundamentais de 
y” +b2y = 0. 


Dependência Linear 


Dizemos que duas funções yı (t) e y2(t) são linearmente dependentes (L.D.) 
em um intervalo I, se uma das funções é um múltiplo escalar da outra, ou 
seja, se 


yı(t) = &y2(t) ou y2(t) =ay(t), paratodot € I. 


Caso contrário, dizemos que elas são linearmente independentes (L.I.). 


Se duas funções são L.D. em um intervalo I, então 


Wim, y2] (t) = det | yl) | =0, paratodot EI, 

2 
pois uma coluna da matriz acima é um múltiplo escalar da outra. Assim, 
vale o seguinte resultado. 


Teorema 4.3. Se yi(t) e ya(t) são funções tais que 


= yı(to) y2(to) 
W[y1, y2] (to) = det | OE #0, paraalgum to € I, 


então y(t) e y2(t) são linearmente independentes (L.I.) em I. Assim, soluções 
fundamentais de (4.1) são L.I. 


Observe que o wronskiano pode ser calculado para quaisquer par de 
funções, mesmo que elas não sejam soluções de uma equação diferencial. 
Também os conceitos de dependência e independência linear são defini- 
dos para duas funções que podem ou não ser soluções de uma equação 
diferencial. 
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Exemplo 4.2. Seja b um número real não nulo. Mostramos no exemplo 
anterior que y1(t) = cos bt e y>(t) = sen bt são soluções L.I. da equação 


y" +by =0. 


A recíproca do Teorema 4.3 não é verdadeira, ou seja, mesmo que 
Wy y2] (t)=0, para todot €R, 


não significa que as funções sejam linearmente dependentes. Vejamos o 
próximo exemplo. 


1 Ty 


0.5 


-T 
af —0.5 0 0.5 1 =} —0.5 0 0.5 1 


Figura 4.1 -n(t) = t? e y2(t) = t|t| são L.I. mas o wronskiano é igual a 
zero para todo t 


2 > 
Exemplo 4.3. Sejam y(t) = P e y2(t) = t|t| = i r - E ; 


Wyry = det[ É fl | =o 
y 2 2t] i 
Apesar de o wronskiano ser zero para todo t € R as funções y1 e y2 são L.I., 


pois uma função não é múltiplo escalar da outra. Para t > 0, y2(t) = y1 (t) 
e para t < 0, y2(t) = —y1 (t). 


4.1.2 Fórmula de Euler 


Queremos definir a função exponencial y(t) = e'! para números complexos 
r = a + ib de forma que satisfaça as propriedades 


elatib)t — enteibt, (4.5) 
é (e”) = re. (4.6) 


Sendo assim, a função z(t) = e™ é solução da equação y” + b?y = 0. Pois 
pela propriedade (4.6) 


z(t) = ibe, z"(t) = —be = —b?z(t) 


e assim 
z” (t) +b?z(t) = 0. 
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Assim, z(t) = e™ é solução do problema de valor inicial 


y" ai b?y = 0, 
y(0) = 1,y'(0) = ib. 


Agora, como mostramos no Exemplo 4.1, que y(t) = cosbt e y2(t) = 
sen bt são soluções fundamentais de y” + b2y = 0, então, pelo Teorema 4.2, 
existem constantes c4 e cy tais que 


z(t) = e! = cycosbt + co senbt. (4.7) 
Vamos determinar estas constantes c4 e c2. Substituindo-se t = O na equação 
(4.7), obtemos que c4 = 1. Derivando a equação (4.7) em relação a t, 
obtemos 

ibe! = —c1b sen bt + cobcosbt. (4.8) 
Substituindo-se t = O na equação (4.8), obtemos que c2 = i. Assim, 


substituindo-se c; = 1 e c2 = i, já obtidos na equação (4.7), obtemos 
et = cosbt + isen bt. 
Portanto, pela propriedade (4.5), 
elit = eatpibt = et (cos bt + isen bt). (4.9) 
Tomando t = 1, temos 
ertib = en( 


cosb +isenb). 


Esta equação é conhecida como fórmula de Euler. 


Exemplo 4.4. Usando a fórmula de Euler temos que 
éT = 1, d? =i, eP?+tii = V2 +iv2, 


que foram obtidas tomando 


7T 7T 
a=0,b=7; a=0,b > a In2, b=, 


respectivamente. 


4.1.3 Obtendo uma Segunda Solução 


Considere uma equação linear de 2º ordem homogênea 


dy 


Ap tay =o. (410) 


dt 

Seja y1 (t) uma solução conhecida da equação acima num intervalo I tal que 
yı(t) £ O para todo t € I. Vamos procurar uma segunda solução da 
equação da forma 


y(t) = v(t)y (t). 
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Como 
vO)=-on+no e voy +2yyo +y", 


então y(t) é solução da equação se, e somente se, 


y” +p(t)y' + a(Oy = y0” +0 (24 + plbya) +o(yi +p(t)yi +4(t)y1) = 0. 


Como y1(t) é solução da equação (4.10), então yY + p(Dy, +4(t)yı = 0, e 
assim a equação anterior se torna 


gro" +v (2y; + p(t)y1) = 0. (4.11) 
Seja w(t) = v' (t). Então a equação (4.11) pode ser escrita como 

yw’ + (2y; + p(t)y1)w =0. 
Esta é uma equação de 1° ordem separável, que pode ser escrita como 


Ee i 


w y1 
Integrando-se, obtemos 
In |w] = -21n |y1| - [o(nát+c, 
que, usando propriedade do logaritmo, pode ser reescrita como 
In jun =—-— / p(t)dt +C. 
Explicitando w(t), obtemos 
em J p(t)dt 


y(t) ` 


Como w(t) = v' (t), resolvendo a equação para v(t): 


em J p(t)dt 
v(t) = af dt + Co. 


w(t) = C4 


y(t)? 
Substituindo-se v(t) em y(t) = v(t)y1 (t), obtemos 
= f plidt 
yO = 00H = Conte) [Coat + Con) 


Tomando-se c2 = 0 e c4 = 1, obtemos uma segunda solução da equação 
(4.10) 
e7 J p(t)dt 


y(t)? i 


Vamos ver que y1 (t) dada e y2(t) obtida por (4.12) são soluções fundamen- 
tais da equação (4.10). 


yl) = (8 | (412) 


Win) = de H pO |- 


- f p(t)dt 
Sd yı(t) nE) Sp dt 
et |, i rp f EOE gpa tia 
nt neS Sear t SO 


= e fPOt zo paratodoteR. 
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Assim, se y1 (t) é uma solução conhecida da equação (4.10) e y2 (t) é dada 
por (4.12), então 

y(t) = c1y1 (t) + coya(t) 
é solução geral da equação (4.10). 


Atenção: Não se deve memorizar a fórmula obtida para y>(t). O que fize- 
mos aqui foi mostrar o caminho que deve ser seguido para encontrar uma 
segunda solução da equação linear homogênea de 2º ordem. 


No próximo exemplo vamos seguir os mesmos passos que seguimos no 
caso geral. 


Exemplo 4.5. Considere a equação 
ay" +by' +cy=0 coma,b,c €R tais que b? — 4ac = Oe a £0. 
Deixamos como exercício verificar que y1 (t) = ez é uma solução da 
equação acima. Vamos procurar uma segunda solução da forma 
—b 
y(t) = v(Dy(t) = v(t)e!,em que r = a” 


Como 
y (t) =v + e yŒ) = v" (tjet +2ro' (tjet +ro(De, 
então y(t) é solução da equação se, e somente se, 


ay" + by' + cy = a(v" +2rv' +r°v) +b(0 — rv) + cv = 


av" + (2ar + b)o' + (ar? +br +c) =0. 
Como r = E é a única solução da equação ar? + br +c = 0 e 2ar +b = 0, 
então 
v” (t) =0. 
Seja w(t) = v' (t). Então a equação v” (t) = O torna-se w'(t) = 0, que tem 
solução w(t) = Cy. Resolvendo a equação v'(t) = w(t) = C1, obtemos 
v(t) =Ct+C e y(t) = (Cit + Co)e'. 


Tomando-se C = 0 e C4 = 1, obtemos 


volt) = te”, 
Vamos ver que yn(t) = e e y2(t) = te", em quer = Ro são soluções 
fundamentais da equação 
nt) val) | et te"! 
ael pa o| T ize ate! 
= 2rt 1 t 
Ra det | r (1+rt) | 


Assim 

y(t) = ce + catet, em quer = E. 
é a solução geral da equação ay” + by' + cy = 0, tal que b? — 4ac = 0 e 
a £40. 
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4.1.4 Equações Homogêneas com Coeficientes Constantes 


Vamos tratar de equações da forma 


d? d 
a b +ey=0, a,bceR, ao. (4.13) 
Para esta equação existem valores constantes de r tais que y(t) = e"! é uma 
solução. 
d d? 
Substituindo-se y = e", = = re" e A = Pe! em (4.13), obtemos 


are + bre"! + cet = (ar? + br +c)e" = 0. 


Como e'! Æ 0, então y(t) = e” é solução de (4.13) se, e somente se, r é 
solução da equação 

ar? +br+c=0 (4.14) 
que é chamada equação característica de (4.13). 


Observe que a equação característica pode ser obtida da equação diferen- 
cial com coeficientes constantes trocando-se y” por r°, y' por r e y por 1. 


Como uma equação de 2° grau pode ter duas raízes reais, somente uma 
raíz real ou duas raízes complexas, usando a equação característica, pode- 
mos chegar a três situações distintas. 


A Equação Característica Tem Duas Raízes Reais 


Se a equação característica de (4.13) tem duas raízes reais (distintas), r4 e 


r2, então 
y(t) =e e y(t) = e 
são soluções fundamentais, pois 
nHO nH] ent e 
der | o) volt) 1 det | ment met 
= entes der | tA | 
ri 12 


= (r —-r)e"it 40, paratodot €R. 


Assim, no caso em que a equação característica tem duas raízes reais dis- 
tintas r1 e ro, 
y(t) = c1" + cze" 


é a solução geral de (4.13). 


Exemplo 4.6. Seja w um número real positivo. Vamos encontrar a solução 
geral da equação 
y" o w2y = 0. 


A equação característica desta equação diferencial é 


"==, 
que tem como raízes r4 = w e r2 = —w. Assim, a solução geral da equação 
diferencial acima é 
y(t) = cre! + ce ™™. 
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Figura 4.2 - Algumas soluções da equação do Exemplo 4.6 


A Equação Característica Tem Somente Uma Raiz Real 


Se A = b? — 4ac = 0, então a equação característica (4.14) tem somente uma 


raiz real r4 = — x Neste caso, 
a 


yilt) = e 


é solução da equação diferencial (4.13). 


No Exemplo 4.5, página 72, mostramos como encontrar uma segunda solução 
para esta equação. Nesse exenplo mostramos que y>(t) = te"! também é 
solução da equação (4.13) e que são soluções fundamentais da equação di- 
ferencial (4.13). 


Portanto, no caso em que a equação característica tem somente uma raíz 


real rı = Er. 


y(t) = cre"! + cote"! 


é a solução geral de (4.13). 


Exemplo 4.7. Vamos encontrar a solução geral da equação 
y” +2y' +y=0. 
A equação característica é 
?+2r+1=0, 
que tem como raiz rı = —1. Assim, a solução geral da equação é 


y(t) = cre! + cate. 
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Figura 4.3 - Algumas soluções da equação do Exemplo 4.7 


A Equação Característica Tem Duas Raízes Complexas 


Se a equação característica de (4.13) tem duas raízes complexas, então elas 
são números complexos conjugados, ou seja, ser; = a + if é uma raiz da 
equação característica (4.14), então a outra raiz é r2 = «a — iB. Neste caso, 
pela fórmula de Euler (4.9): 


y(t) = et = e®tiP = et (cos Bt +isenpt) e 
yalt) = et = eib) = et (cos(—Bt) + isen(—Bt)) = e“ (cos Bt — isen Bt). 


Pela análise feita no início desta seção, sabemos que y1 (t) e y>(t) são soluções 
complexas da equação diferencial (4.13). Além disso, 


yi(t) 
iii | vil) 


=== 
II 


det | 
ente der | 1 1 | 
ri 


= (r-r)etitot = (—2iß)e™ +0, para todot €R, 


N~ 

A 
+ o —+ 
so n— 


ou seja, yi (t) e y>(t) são soluções fundamentais de (4.13). Assim, no caso 
em que a equação característica tem duas raízes complexas rı = a + if e 
T2 = & — if, 

y(t) = Ce! +C, CC EC 
é a solução geral complexa de (4.13). 


Vamos encontrar um conjunto fundamental de soluções reais. A solução 
geral complexa pode ser escrita como 


y(t) = Cyele+ib)t di Czel*—iP)t 
= Ce" (cos Bt + isen Bt) + Coe" (cos Bt — i sen Bt) 
= (C1 + C2)e™ cos Bt + i(C1 — C2)e“ sen Bt. (4.15) 


1 
Tomando C4 = € = 5 em (4.15), temos a solução real 


u(t) = e" cosßt. 
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1 1 
Tomando C4 = > eC = — z temos a solução real 


v(t) = e“ sen pft. 


Vamos mostrar, agora, que, se as raízes da equação característica são com- 
plexas, então u(t) e v(t) são soluções fundamentais de (4.13). 


a | E. | 


det e"! cos Bt e“! sen Bt 
j e*t (acos pt — psen Bt) e% (œ sen Bt + Bcos Bt) 


ou (uam [ SER nB | egaa] ce, nA |) 


peZrt * 0, paratodot €R. 


Assim, no caso em que a equação característica tem duas raízes complexas, 
m=a+tiber,=a-iB, 


y(t) = cet cos Bt + coe"! sen Bt 
é a solução geral de (4.13). 
Exemplo 4.8. Seja w um número real positivo. Vamos encontrar a solução 
geral da equação 
y! +w?y =0. 


A equação característica desta equação diferencial é 


+w? =0, 
que tem como raízes rm = iw e rọ = —iw. Assim, a solução geral da 
equação diferencial acima é 
y(t) = cy cos wt + cp sen wt. (4.16) 


Escrevendo o par (c1, c2) em coordenadas polares, temos que 


(cy, C2) 


Rcosó, 


Rsenô. (417) 


» 
Ve 
C3r «CS 
N e 
|l 


Substituindo-se os valores de cy e c2 na equação (4.16), obtemos 


y(t) = R (cos ô cos (wt) + sen ô sen (wt)) = Rcos(wt — ô), 


em que R = ga + c2 e ô são obtidos de (4.17). 
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y(t) = Rcos(wt- 5),8>0,0>0, R>0 


Figura 4.4 - Uma solução da equação do Exemplo 4.8 
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. (a) Mostre que y (x) = x? e y2 (x) =x 


Exercícios 


. Mostre que y1(x) = x? é solução da equação diferencial 


2x°y" = xy —9y=0, 


Encontre uma função u(x) tal que y>(x) = u(x)y (x) seja solução da equação dada. Prove que 
as duas soluções yı (x) e y>(x) são soluções fundamentais. 


. Mostre que y1 (x) = x71, x > 0, é solução da equação diferencial 


e Bay du O. 


Encontre uma função u(x) tal que y>(x) = u(x)y (x) seja solução da equação dada. Prove que 
as duas soluções yı (x) e y>(x) são soluções fundamentais. 


. Encontre a solução geral da equação 


y” +2y' + ay = 0, 


para a > 1, para « = 1 e para 4 < 1. 
3 


. (a) Determine qual ou quais das funções zı (x) = x?, z2(x) = x? e z3(x) = e * são soluções da 


equação 


(x +3)y" +(x +2)y' -y =0. 
(b) Seja y (x) uma das soluções obtidas no item anterior. Determine uma segunda solução 
y2(x) de forma que y; (x) e y>(x) sejam soluções fundamentais da equação. 
(c) Determine a solução geral da equação 
(x +3)y” + (x+2)y' -y =0 
e obtenha a solução do problema de valor inicial 


| (x +3)y” +(x +2)y' -y =0, 


y) =1, 
y'(1)=3. 
Justifique sua resposta! 
3 são soluções da equação 
x?y" = 6xy' + 10y = 0, 
(b) Obtenha a solução do problema de valor inicial 
x2y" — 6xy' + 10y = 0, 
"O 
y (1) =3. 


Justifique sua resposta! 


. Mostre que a solução do problema y” + 2y' = 0,y(0) = a, y' (0) = b tende para uma constante 


quando tł — +00. Determine esta constante. 


. Mostre que, se 0 < b < 2, então toda solução de y” + by’ + y = 0 tende a zero quando t — +00. 


. Considere o problema y” — 4y = 0,y(0) = 0,y'(0) = b # 0. Mostre que y(t) £ O para todo 


t0. 
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9! 


10. 


ML, 


12 


13. 


14. 


15. 


Considere o problema y” — y' + y = 0, y(0) = 2,y' (0) = b. Determine os valores de b para os 
quais a solução y(t) — +% quando t — +00. 


Considere a equação y” + 2by' +y = 0. Para quais valores de b a solução y(t) tende a zero 
quando t — +0, independentemente das condições iniciais. 


As equações de Euler são equações que podem ser escritas na forma 


2 
pa Y + bx L Fey = 0, (4.18) 


em que b e c são constantes reais. 


Mostre que existem valores constantes de r tais que y(x) = x” é uma solução de (4.18). Além 
disso, mostre que y(x) = x” é solução da equação (4.18) se, e somente se, 


P+(b-1)r+c=0. (4.19) 
A equação (4.19) é chamada equação indicial de (4.18). 
Mostre que se a equação indicial (4.19) tem duas raízes reais (distintas), rj e r2, então 
y(x) =x" e ya(x) =x? 
são soluções fundamentais de (4.18) e, portanto, 
ylx) = cx + cox? 


é a solução geral de (4.18), para x > 0. 


Se a equação indicial (4.19) tem duas raízes complexas, rı = a + iß e r2 = a — ip, use a fórmula de 
Euler para escrever a solução geral complexa em termos das soluções reais, para x > 0, 


u(x) = x*cos(Blnx) e v(x) =x"sen(Blnx). 
Mostre que estas soluções são soluções fundamentais de (4.18) e, portanto, 
y(x) = c1x" cos( ln x) + c2x“ sen(BIn x) 
é a solução geral de (4.18), para x > 0. 


1—b 
Se a equação indicial (4.19) tem somente uma raíz real, m = F determine uma segunda 


solução linearmente independente da forma y2 (x) = v(x)yı (x) = v(x)x T, para x > 0. Mostre 
1-b 1=b : 

que y1(x) = x 7 ey>(x) =x Inx são soluções fundamentais de (4.18) e, portanto, a solução 

geral de (4.18), para x > 0, é 


1-b 1-b 
y(x) =cıx Z? +cx7 lnx. 


Use os exercícios anteriores para encontrar a solução geral das seguintes equações: 


(a) xy" +4xy' + 2y = 0. 
(b) x?y" — 3xy' + 4y = 0. 
(c) x2y" +3xy' + 5y = 0. 


TO 


Equações Diferenciais Lineares 
de 2º Ordem 


2° Parte 


5.1 Equações Não Homogêneas 


OBJETIVOS 


Ao terminar esta seção você deverá ser capaz de: 


e Identificar uma equação diferencial linear de 2º ordem não homogênea. 
e Encontrar a solução geral de uma equação diferencial linear de 2º ordem com 
coeficientes constantes não homogênea. 


e Resolver um problema de valor inicial correspondente a uma equação diferencial 
linear de 2? ordem não homogênea. 


Uma equação diferencial linear de 2? ordem não homogênea é uma equação 
da forma 


(GE + (E + el = 0, 


com f(t) uma função não nula, que pode ser reescrita como 


dy 

t TO gr T Cy = FC). (5.1) 
Teorema 5.1. Seja yp(t) uma solução particular da equação (5.1). Sejam yı (t) 
e y2(t) soluções fundamentais da equação homogênea correspondente. Então a 
solução geral da equação não homogênea (5.1) é 


y(t) = yp(t) + c1y1 (t) + coyo(t). 


Demonstração. Seja y(t) uma solução qualquer de (5.1) e yp(t) uma solução 
particular de (5.1), então Y(t) = y(t) — yp(t) é solução da equação ho- 
mogênea associada 

dy dy = 5.2 
a trO tI =O (5.2) 
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Pois, 


Assim, se y(t) e y>(t) são soluções fundamentais da equação homogênea 
associada (5.2), existem constantes cj e co tais que 


Y(t) = y(t) — yp(t) = c1y1 (t) + coyolt), 


ou seja, se y(t) é uma solução qualquer de (5.1) e y1 (t) e y2(t) são soluções 
fundamentais da equação homogênea associada (5.2), então 


y(t) = c1y1 (t) + c2y2 (t) + yp (t). (5.3) 
E 


Portanto, para encontrar a solução geral de uma equação linear de 2º or- 
dem não homogênea precisamos encontrar uma solução particular e duas 
soluções fundamentais da equação homogênea correspondente. 


Teorema 5.2 (Princípio da Superposição para Equações Não Homogêneas). 
Se yo? (t (t) é uma solução de 


y” + p(Dy + q(Dy = fi (t) 


e y® (t) é uma solução de 


y” + p(y + g(Dy = fat), 
então, yp(t) = yP (t) + yP (t) é solução de 
y” + p(Dy + a(y = filt) + fab). 
Demonstração. 


yp)" + py, pO alyth = 
= WP (DHP D)" + POD Hyp OY HOU O Hyp (1) = 
= yp (H+ tO nus (H Hyp ("+ pO (D +a (A) = 

=f(t) =fo(t) 


= f(t)+fott), 
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pois y (t) é solução da equação 


y" + py +4(t)y = fi (t) 
e y” (t), da equação 


y” +p(t)y' +q(t)y = folt). 


Exemplo 5.1. A função yı (t) = ; é solução da equação diferencial 
y” +4y=t 
(verifique!) e a função y>(t) = : sen(2t) é solução da equação 
y” + 4y = 2 cos(2t) 


(verifique!). Pelo Princípio da Superposição para Equações Não Homogêneas 


t t 
(Teorema 5.2) y(t) = Z + > sen(2t) é solução da equação 


y” + 4y = 2 cos(2t) + t. 
5.1.1 Equações Não Homogêneas com 
Coeficientes Constantes 
Vamos tratar de equações da forma 


dy dy 
a2 F st), 


em que a, b e c são números reais, a 0. 


Método dos Coeficientes a Determinar 


Este método funciona quando a função g(t) tem uma das seguintes formas: 


(1) g(t) = ao +... + ant”, em que aọ,...,an ER. 
Neste caso deve-se procurar uma solução particular da forma 


yp(t) = È (Ao +... + Ant”), 


em que s é o menor inteiro não negativo que garanta que nenhuma 
parcela de yp(t) seja solução da equação homogênea correspondente, 
e Ag, .. . , An são coeficientes a ser determinados substituindo-se yp(t) 


na equação (5.4). O Exemplo 5.2 ilustra este caso. 


(2) g(t) = (ag +... + ant” )e™, em que ag,...,ap,a ER. 
Neste caso deve-se procurar uma solução particular da forma 


volt) = Ë (Ao +... + Ante, 


em que s é o menor inteiro não negativo que garanta que nenhuma 
parcela de yp(t) seja solução da equação homogênea correspondente 
e Ag, .. -, An são coeficientes a ser determinados substituindo-se yp (t) 


na equação (5.4). O Exemplo 5.3 ilustra este caso. 
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(3) f(t) = (ao +... + ant” je" cos Bt ou f(t) = (ao +... + ant" Je" sen Bt, 
em que 4o, ..., an, &, b ER. 
Neste caso deve-se procurar uma solução particular da forma 


yp(t) = EI(Ao +... + Ant” )e" cos Bt + (Bo +... + Bat" )e sen Bt), 


em que s é o menor inteiro não negativo que garanta que nenhuma 
parcela de yp(t) seja solução da equação homogênea correspondente 
e Ag,..., An, Bo, . - ., Bn são coeficientes a serem determinados subs- 
tituindo-se yp(t) na equação (5.4). O Exemplo 5.4 ilustra este caso. 


Exemplo 5.2. Vamos encontrar a solução geral da equação 
y! +y =2+P. 


Precisamos encontrar a solução geral da equação homogênea correspon- 
dente y” +y’ = 0. A equação característica é 


P+r=0, 


que tem como raízes r4 = 0 e r2 = —1. Assim, a solução geral da equação 
homogênea correspondente y” +y’ = 0 é 
y(t) = c1 + ce. 


O segundo membro da equação g(t) = 2 + £ é da forma (1). Vamos pro- 
curar uma solução particular da forma 


yp(t) = t (Ao + A1t + Aat?) = Aot + Ag + A2? 


O valor de s é igual a 1, pois para s = 0, a parcela Ag é solução da equação 
homogênea (c2 = 0 e cy = Ag). 


Volt) = Ao +2A1t +342? 
yp(t) = 244 + 6At. 
Substituindo y, (t) e y, (t) na equação y” +y’ = 2 + t2, obtemos 


(2A1 +64A2t) + (Ao + 2A44t +342t) = 


(Ao +241) + (241 + 64A2)t +342 =2 +t 


Comparando os termos de mesmo grau obtemos o sistema linear 


Ao + 2A = 2 
2A; + 64 = 0 
34, = 1, 
que tem solução Ao = 4, Aı = —1 e 4 = 1/3. Assim, uma solução 


particular da equação não homogênea é 
Yp(t) -g4-P4lg 
Í 3 


e a solução geral da equação não homogênea é 


1 
y(t) = c1 +ee™ + 4t -t +P. 
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Figura 5.1 - Algumas soluções da equação do Exemplo 5.4 


Exemplo 5.3. Vamos encontrar a solução geral da equação 
y” +2y' +y = (2+). 


Precisamos encontrar a solução geral da equação homogênea correspon- 
dente y” +2y' +y = 0. A equação característica é 


r +2r+1 =0, 


que tem como raiz r4 = —1. Assim, a solução geral da equação homogênea 
correspondente y” +2y' +y =0é 


y(t) = cre! + czte™. 


O segundo membro da equação g(t) = (2 + t)e™* é da forma (2). Vamos 
procurar uma solução particular da forma 


vp(t) = P (Ao + A1t)e™ = (Ag + Ate", 


O valor de s é igual a 2, pois para s = 0 as parcelas Age”! e Aste | são 
soluções da equação homogênea (cı = 49, c2 = 0 e cı = 0, = As)e 
para s = 1 a parcela Apte™* é solução da equação homogênea (cı = 0 e 


Cc = Ao). 


y(t) = (240t H (344 — Ag)f? At) et 


y(t) = (240 + (641 — 4A0)t + (Ao 64) +41) et, 


[4 
j 


y = (2 + t)e™t, obtemos 


Substituindo y, (t) e yp(t) na equação y” + 2y 


(240 + (641 — 440)t + (Ao — 6A1)Ë A At) e~t + 


$2 (240t (341 — Ao)t? At) et + 
+ (Ao Ando = Qro 
Simplificando o primeiro membro, obtemos 
(249 +64t)e !=(24+te! > 249+64t=2+4t. 
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Comparando os termos de mesmo grau, obtemos o sistema linear 


249 = 2 
6A1 = 1, 
que tem solução Ao = 1 e 4; = 1/6. Assim, uma solução particular da 


equação não homogênea é 


1 £ 
volt) = (P+ Ge! 


e a solução geral da equação não homogênea é 


y(t) = qe! + czte™t + (P + Pe. 


Figura 5.2 - Algumas soluções da equação do Exemplo 5.2 


Exemplo 5.4. Vamos encontrar a solução geral da equação 
y” +2y' + 2y = é cost. 


Precisamos encontrar a solução geral da equação homogênea correspon- 
dente y” + 2y' + 2y = 0. A equação característica é 


r? +2r+2=0, 


que tem como raízes rı = —1 +i er? = —1 — i. Assim, a solução geral da 
equação homogênea correspondente y” + 2y' + 2y = 0 é 


y(t) = c1e™* cost + cze™* sent. 


O segundo membro da equação g(t) = e' cost é da forma (3). Vamos pro- 
curar uma solução particular da forma 


yp(t) = tº(Ae' cost + Be' sent) = Ae' cost + Be' sent. 


O valor de s é igual a 0, pois nenhuma parcela de yp(t) é solução da 
equação homogênea. 


yp(t) = A(e' cost — e' sent) + B(e'sent + e cost+) = 
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(A+ B)e' cost + (B — A)e sent 
vp(t) = 2Be' cost — 2Ae' sent. 
Substituindo yp(t) e y, (t) na equação y” + 2y' +y = e cost, obtemos 
2Be' cost — 2Ae sen t+ 
+2((A+B)e'cost+ (B — A) sent) + 
+2(Ae'cost+Be'sent) = e cost. 


Simplificando o primeiro membro, obtemos 
(4A +4B)e' cost + (4B — 4A)e sent = é cost. 


Comparando os coeficientes de e' cost e de ef sent, obtemos o sistema li- 
near 
4A + 4B = 1 
i —4A + 4B = 0, 


que tem solução A = 1/8 e B = 1/8. Assim, uma solução particular da 
equação não homogênea é 


(t) = gaia seg 
Yp “8 8 


e a solução geral da equação não homogênea é 


t 


1 
y(t) = c1e™* cos t + coe ! sent + ze (cos t+sent). 


Figura 5.3 - Algumas soluções da equação do Exemplo 5.4 
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1. Encontre a solução geral das equações: 
(a) y +5y + 6y = xe ™™., 
(b) y” — 4y' + 6y = 3x. 
(c) y” +4y = Aai 
(d) y” +2y = # +2. 


2. Resolva os problemas de valor inicial: 
(a) y” +y'-2y=?+3, y(0)=0, y'(0)=0. 
(b) r +2y' +y =3sen(2t), y(0)=0, y(0)=0. 
(c) y” —4y' +4y=3e™*, y(0)=0, y'(0)=0. 
(d) 2y" +2y'+y =, y(0)=0, y'(0)=0. 


3. (a) Encontre a solução geral da equação 


y” +2y' +ay=0, 


para « > 1, para « = 1 e para «& < 1. 
(b) Determine a forma adequada para uma solução particular da equação 


y” +2y' + ay = te“ sen(vVa — 1t), 


para «a > 1. 
(c) Para quais valores de « todas as soluções tendem a zero quando t — -+oo. 


Equações Diferenciais Lineares 
de 2º Ordem 


3º Parte 


6.1 Oscilações 


OBJETIVOS: 
Ao terminar esta seção você deverá ser capaz de: 


e Modelar um problema envolvendo um sistema massa-mola. 
e Resolver um problema de valor inicial correspondente ao modelo encontrado. 


e Identificar o tipo de movimento: oscilatório (periódico, quase-periódico) e não 
oscilatório. 


e Calcular, no caso de movimento periódico, o período, a frequência e a fase. 


e Calcular, no caso de movimento quase-periódico, o quase-período, a quase- 
-frequência e a fase. 


e Identificar as condições para ressonância e saber calcular a frequência de 
ressonância. 


e Calcular a solução estacionária e a solução transiente, no caso de oscilações 
forçadas amortecidas. 


Considere um sistema massa-mola na vertical. Seja L o alongamento pro- 
vocado na mola pela colocação da massa m quando o sistema está em 
equilíbrio. Neste caso a magnitude da força elástica é igual à magnitude 
do peso, ou seja, 

mg = kL. (6.1) 


Aqui k é chamada constante da mola. Seja y(t) o alongamento da mola em 
um instante t. Defina a nova função 


u(t) = y(t) - L. 
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Sobre a massa agem o seu peso, 
P = mg, 


a força da mola que é proporcional ao seu alongamento e tem sentido 


oposto a ele, 
Fe = —ky(t) = —k(u(t) + L), 


uma força de resistência proporcional à velocidade, 
F, = —yy' (t) = — yu’ (t) 
e uma força externa Fex+. Aqui y é a constante de amortecimento. 


Pela segunda lei de Newton, temos que 


my" (t) = mg — ky(t) — yy’ (t) + Fext 


ou escrevendo em termos de u(t) = y(t) — L: 


mu” (t) = mg — K(L + u(t)) — yu’ (t) + Eext. (6.2) 
Assim, por (6.1) e (6.2), u(t) satisfaz a seguinte equação diferencial 
mu” (t) + qyu'(t) + ku(t) = Fext. (6.3) 


que é a mesma equação que satisfaz x(t), no caso de a mola estar na posição 
horizontal. Verifique! 


E 


Figura 6.1 - Sistema Massa-mola na vertical 
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6.1.1 Oscilações Livres 
Sem Amortecimento 


Como as oscilações são livres, Fex, = 0 e como são não amortecidas, y = 0. 
Assim, a equação (6.3) para o movimento da massa é 


mu” + ku = 0. 


A equação característica é 


Jk 
mr +k=0 < r=g i. 
m 


Assim, a solução geral da equação é 


u(t) = cy cos CE ) + c2 sen (v ) ; 


Seja wo = k, Então a equação acima pode ser escrita em termos de wo 
como 
u(t) = cy cos (wot) + c2 sen (wot) . (6.4) 
Marcando o ponto (c1, c2) no plano e escrevendo em coordenadas polares 
temos que 
y 


(C1, C2) 


C2 


o = Rcosó, (65) 


co = Rsenó. 


Substituindo-se os valores de cy e c2 obtidos de (6.5) na equação (6.4), 
obtemos 


u(t) = Rcosôcos (wot) + R sen ôsen (wot) 
R (cos ô cos (wot) + sen ô sen (wot) ) 
= Rcos(wot-— ô), 


Aqui foi usada a relação 
cos(a — b) = cosa cos b + sen asen b. 
wọ é chamada frequência natural do sistema, ô a fase e R a amplitude. 


sda 27 

Neste caso a solução da equação é periódica de período T = D Este 
0 
movimento oscilatório é chamado movimento harmônico simples. 
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xT 


Figura 6.2 - Sistema massa-mola livre não amortecido 


hu ult) = R osloni = 5 


wd A ye — ô) 
w=\ 5 


Figura 6.3 - Solução do sistema massa-mola livre não amortecido 


Exemplo 6.1. Sabendo-se que o problema de valor inicial que descreve um 
sistema massa-mola é dado por 


y” +5y=0, y(0)=1, y(0)=0. 


(a) Encontre a solução geral da equação diferencial e resolva o problema 
de valor inicial. Determine a amplitude, a frequência, a fase e o 
período. 


(b) Esboce o gráfico da solução obtida. 
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Solução: 
(a) A equação característica é 


r? +5=0, 


que tem como raízes r = +/5i. Assim, a solução geral da equação é 
y(t) = cy cos (v5 t) + c2 sen (v5 t) i 


Para resolver o problema de valor inicial precisamos calcular a deri- 
vada da solução geral 


y' (t) = -v5 c1 sen (51) + v5 c2 cos (51) . 


Substituindo-set = 0, y = 1 e y' = 0, obtemos cı = 1e c = 0ea 
solução do problema de valor inicial é 


y(t) = cos (v5 t) ; 


A amplitude é igual a 1, a frequência é igual a V5, a fase é igual a 
zero e o período é igual a 27r / v5. 


dy 


Com Amortecimento 
Como as oscilações são livres, Ex, = 0. Assim, a equação (6.3) para o 
movimento da massa é 
mu” + yu’ + ku =0 
A equação característica é mr? + yr +k = 0 e A = 9? — 4km 
Aqui temos três casos a considerar: 


(a) Se A = 9°? — 4km > 0 ou y > 2v'km, neste caso 


u(t) = cet + cze", 


em que 
= =E A = L y2 — 4km 
no- EVA vt vY "0, 
' 2m 2m 
Este caso é chamado superamortecimento e a solução u(t) — 0 quando 
t — +0. 
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(b) Se A = 9º — 4km = 0 ou y = 2Vkm, neste caso 


zi = 
u(t) = cje 2m + cote” m 


Este caso é chamado amortecimento crítico e a solução u(t) > 0 
quando t — “oo. 


(c) Se A = 9? — 4km < 0 ou 0 < y <2wkm, neste caso 


u(t) = eim (cy cos ut + c2 sen yt), (6.6) 


em que 


U= ima A0: 


a 2 ia 


2 
Aqui, 4 é chamado quase-frequência e T = ZT é chamado quase- 


-período. Escrevendo novamente o par (c1, c2) em coordenadas po- 
lares temos que 


(c1, C2) 


C2 


Rcosó, 


Co o = 
j E = Rsenô. e) 


Substituindo-se os valores de cj e c} na equação (6.6), obtemos 


u(t) = em (R cos ô cos (ut) + R sen ffi sen (ut)) = Rem cos(ut — ô), 


em que R = Vê + E e ô são obtidos de (6.7). 
Este caso é chamado subamortecimento e a solução u(t) — 0 quando 


p ; Said da ; x 
t — +o. Este é um movimento oscilatório com amplitude Re” 2m, 
chamado quase-periódico. 


Observe que nos três casos a solução tende a zero quando t tende a +00. 


AULA 6 


= 
0 x 


Figura 6.4 - Sistema massa-mola livre com amortecimento 


u u(t) = elit + ce'2! 
-y + Vy? — 4km 
[1,2 = = A 
2m 


C + C2 = Uo 


Figura 6.5 - Algumas soluções do sistema massa-mola livre com supera- 


mortecimento 
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Figura 6.6 - Algumas soluções do sistema massa-mola livre com amorteci- 
mento crítico 


Sub Amortecimento 


yt 
2m (C14 COS ut + C2 sen pt) 


Figura 6.7 - Algumas soluções do sistema massa-mola livre com subamor- 
tecimento 
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Figura 6.8 - Solução típica do sistema massa-mola livre com subamorteci- 
mento 


super amortecimento, y > 2 Ykm 


/ 


amortecimento crítico, y = 2 vkm 


sub amortecimento, y < 2 vkm 


Figura 6.9 - Comparação das soluções do sistema massa-mola livre com 
amortecimento para diferentes valores da constante de amortecimento y 
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6.1.2 Oscilações Forçadas 


Vamos supor que uma força externa periódica da forma Fx, = Focos(wt) 
seja aplicada à massa. Então a equação (6.3) para o movimento da massa é 


mu” + yu! + ku = Focos(wt). 


Oscilações Forçadas sem Amortecimento 


Neste caso a equação diferencial para o movimento da massa é 


mu” + ku = Fycos(wt). (6.8) 


Sabemos que as soluções são da forma 
u(t) = c1 cos (wot) + co sen (wot) + up(t), 
em que, pelo método das constantes a determinar, 
up(t) = t [A cos(wt) + Bsen(wt)] 


é uma solução particular e s é o menor inteiro não negativo que garanta 
que nenhuma parcela de up(t) seja solução da equação homogênea corres- 
pondente e A e B são coeficientes a serem determinados substituindo-se 
up(t) na equação diferencial (6.8). 


Temos dois casos a considerar: 


(a) Se w £ wo. Neste caso s = 0, pois nenhuma das parcelas de up(t) 
é solução da equação homogênea correspondente. Então a solução 
particular é da forma 


up(t) = Acos(wt) + Bsen(wt) 
e a solução geral da equação é da forma 
u(t) = cy cos (wot) + c2 sen (wot) + A cos(wt) + Bsen(wt). 


Deixamos como exercício para o leitor verificar que substituindo-se 
up(t) na equação diferencial (6.8) encontramos 


Fo 


A=; 
m(w — w?) 


e B=0. 


Assim, 


Fo 


Sa) cos(wt). 


u(t) = cy cos (wot) + c2 sen (wot) + 


Neste caso a solução u(t) é oscilatória e limitada. 


(b) Se w = wọ. Neste caso s = 1, pois para s = 0 as parcelas, A cos(wot) 
e Bsen(wot), de up(t), são soluções da equação homogênea corres- 
pondente. Então, a solução particular é da forma 


up(t) = [A cos(wt) + Bsen(wt)] 
e a solução geral da equação é da forma 


u(t) = cy cos (wot) + c2 sen (wot) + t[A cos(wot) + Bsen(wot)). 
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Deixamos como exercício para o leitor verificar que, substituindo-se 
up(t) na equação diferencial (6.8), encontramos 


Ay é Best. 
2mwo 
Assim, 
EEE E EE E) 
= cy cos (wo > sen (wo di sen(wot). 


Neste caso u(t) é oscilatória, mas fica ilimitada quando t tende a +00. 
Este fenômeno é conhecido como ressonância e a frequência w = wo 
é chamada frequência de ressonância. 


Fê = F cos(œt) 


===" 
F =-kx 
MAAA E 


T Lagi 
XxX 


Figura 6.10 - Sistema massa-mola forçado sem amortecimento 


Exemplo 6.2. Vamos considerar o problema de valor inicial 
mu” + ku = Fycos(wt), 
u(0) = 0,u'(0) = 0. 


Temos dois casos a considerar: 


(a) Se w Æ wo. A solução geral da equação é 


Fo 


u(t) =c wot C wot) + ——s 
(t) = c1 cos (wot) + c2 sen (wot) mau) 


cos(wt). 


Derivando e substituindo-se t = 0, u = 0 e u” = 0, obtemos que 


Fo 
m(w3 — w?) 


c1 = co = 0. 


Assim, a solução do problema de valor inicial é 


ij ==" 


VER (cos(wt) — cos(wot)) . 


Como 
cos(A — B) — cos(A + B) = 2sen A sen B 
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então 2F 
u(t) = —— + sen (w1t) sen(wst) 
m(wo — w?) 
em que 
Wo — w w wo +w 
w = = 
1 2 , 2 2 

Como w; = “2 é menor do que w = vos “ então o movimento 


é uma oscilação de frequência w2 com uma amplitude também osci- 


Sn dá 2H 
latória R(t) = E 


é chamado batimento. 


sen(w1t) de frequência w1. Este movimento 


(b) Se w = wo. A solução geral da equação é 


Fo 
2mw 


u(t) = cy cos (wot) + c2 sen (wot) + tsen(wot). 
0 


Já vimos que neste caso u(t) fica ilimitada quando t tende a +œ, que 
é o fenômeno da ressonância. Derivando e substituindo-se t = 0, 
= 0e u’ = 0, obtemos que 


c =0, c0=0 


Assim, a solução do problema de valor inicial é 


F 
u(t) = a tsen(wot). 


Este movimento é uma oscilação de frequência wo com uma ampli- 


tude R(t) = an t que aumenta proporcionalmente a t. 
u u(t) = R sen(wt) sen(w2t), 
R= =i 
m (wZ— w2)” 
w = 2, w = Zoi 


DEPEN 


R senlo 9 > “4 i 
Vs 


Figura 6.11 - Solução do sistema massa-mola, para u(0) = u'(0) = 0, no 
caso de batimento 
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Figura 6.12 - Solução do sistema massa-mola, para u(0) = u'(0) = 0, no 
caso de ressonância 


Oscilações Forçadas com Amortecimento 


mu” + yu! + ku = Focos(wt) 
Seja u(t) = cyu(t) + cou>(t) a solução da equação homogênea correspon- 
dente. Então a solução geral desta equação é 


u(t) = cyu (t) +couo(t) + up(t), 


em que up(t) é uma solução particular. Pelo método das constantes a 
determinar 
up(t) = Ag cos(wt) + Bo sen(wt). 


Substituindo-se up(t) e sua derivada na equação, encontramos 


o Fom (wa — w?) By = Foyw 


A Fá E 
) A A 


em que A = m?(w5 — w)? + y2w?. Podemos escrever 


up(t) = Agcos(wt) + Bosen(wt) = Rcos(wt — ô), 


em que R = 4/Aĝ + Bĝ e ô é tal que Ao = Rcosô e By = Rsenô. Neste 


caso E 
mw — w 
R = =, ô= arccos ag ET 


a 
Vs Ja 


Assim, a solução geral da equação é 
u(t) = c1u1 (t) + c2u2 (t) + Rcos(wt — ô). 


A solução geral da equação homogênea correspondente, c1u1 (t) + c2u2(t), 
é a solução do problema de oscilação livre amortecida e já mostramos que 
tende a zero quando t tende a +0, por isso é chamada solução transiente, 
enquanto a solução particular, R cos(wt — ô), permanece, e por isso é cha- 
mada solução estacionária. 


u(t) = ciu (t) + cous(t) + Rcos(wt — ô) = Rcos(wt — ô), 


para t suficientemente grande. 
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F; =-yv F = F cos(ot) 
F =-kx 
e 


„Fox = F cos(ot) 


F =-yv F E F cos(ot) 


exi 


= 
X 


Figura 6.13 - Sistema massa-mola forçado com amortecimento 


Figura 6.14 - Solução do sistema massa-mola forçado com amortecimento 
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Exercícios 


. Sabendo-se que o problema de valor inicial que descreve um sistema massa-mola é dado por 
y'+5y=0, y(0)=1, y'(0)=0 


(a) Encontre a solução geral da equação diferencial e resolva o problema de valor inicial. De- 
termine a amplitude, a frequência, a fase e o período. 


(b) Esboce o gráfico da solução obtida. 


. Sabendo-se que o problema de valor inicial que descreve um sistema massa-mola é dado por 
2y" +3y=0, y(0)=1, y(0)=0 


(a) Encontre a solução geral da equação e resolva o problema de valor inicial. Determine a 
amplitude, a frequência, a fase e o período. 


(b) Esboce o gráfico da solução obtida. 


. Uma mola, de um sistema massa-mola sem amortecimento, tem constante de elasticidade igual 
a3 N/m. Pendura-se na mola uma massa de 2 kg e o sistema sofre a ação de uma força externa 
de 3 cos(3t). Determine a função que descreve o movimento da massa em qualquer instante t, 
considerando a posição inicial igual ug e a velocidade inicial ug. 


. Se um sistema massa-mola com uma massa de 2 kg e uma mola com constante de elasticidade 
igual 0,5 N/m é colocado em movimento, no instante t = 0, num meio em que a constante 
de amortecimento é igual a 1 N.s/m, determine a posição da massa em qualquer instante t, 
considerando a posição inicial igual a ug e a velocidade inicial uç. 


. Uma massa de 100 gramas estica uma mola 10 centímetros. Suponha que não haja amorteci- 
mento e que a aceleração da gravidade seja de 10º centímetros por segundo ao quadrado. En- 
contre a frequência, o período e a amplitude do movimento. Determine a posição u em função 
do tempo t e faça um esboço do seu gráfico: 


(a) Se a massa é colocada em movimento a partir da sua posição de equilíbrio com uma veloci- 
dade apontada para cima de 4 centímetros por segundo. 

(b) Se a massa é puxada para baixo contraindo a mola 1 centímetro e depois colocada em mo- 
vimento com uma velocidade para baixo de 10 centímetros por segundo. 


(c) Se a massa é puxada para baixo contraindo a mola 2 centímetros e depois é solta. 


. Uma massa de 100 gramas estica uma mola 10 centímetros. A massa está presa a um amorte- 
cedor viscoso. Suponha que a aceleração da gravidade seja de 10° centímetros por segundo ao 
quadrado. 


(a) Para quais valores da constante de amortecimento y o sistema é superamortecido, tem um 
amortecimento crítico e é subamortecido. 

(b) Suponha que o amortecedor exerça uma força de 10º dinas (=gramas-centímetros por 
segundos?) quando a velocidade é de 10 centímetros por segundo. Se a massa é puxada 
para baixo 2 centímetros e depois é solta, determine a posição u em função do tempo t e 
faça um esboço do seu gráfico. Qual o valor do quase período? 
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Uma massa de 100 gramas estica uma mola 10 centímetros. Suponha que não haja amorteci- 
mento e que a aceleração da gravidade seja de 10? centímetros por segundo ao quadrado. Se 
o sistema é colocado em movimento com uma força externa de 9600 cos(6t) dinas, determine a 
posição da massa como função do tempo e faça um esboço do seu gráfico. 


. Uma massa de 100 gramas estica uma mola 10 centímetros. Suponha que não haja amorte- 


cimento e que a aceleração da gravidade seja de 10? centímetros por segundo ao quadrado. 
Se o sistema é colocado em movimento na posição de equilíbrio com uma força externa de 
1000 cos(wt) dinas, para w igual à frequência de ressonância, determine a posição da massa 
como função do tempo e faça um esboço do seu gráfico. 


. Uma massa de 100 gramas estica uma mola 10 centímetros. A massa está presa a um amorte- 


cedor viscoso. Suponha que a aceleração da gravidade seja de 10º centímetros por segundo ao 
quadrado. Suponha que o amortecedor exerça uma força de 4200 dinas quando a velocidade é 
de 1 centímetro por segundo. Se a massa está sob a ação de uma força externa de 26000 cos(6t) 
dinas, determine a posição u em função do tempo t e faça um esboço do seu gráfico, conside- 
rando somente a solução estacionária. 


Considere a equação diferencial do sistema massa-mola forçado sem amortecimento 
mu” + ku = Fycos(wt). 


Mostre que a solução geral: 
(a) Se w £ wo é dada por 


Fo 
E) = t t t). 
u(t) = cy cos (wot) + c2 sen (wot) + e cos(wt) 
(b) Se w = wo é dada por 
F 
u(t) = cy cos (wot) + c2 sen (wot) + mT tsen(wot). 
0 


Mostre que a solução do PVI 


mu” + ku = Fo cos(wt), 
u(0) = 0, u’ (0) = 0. 


(a) Se w £ wo é dada por 


“O a a (e ent) 
(b) Se w = wọ é dada por 
fi a 
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